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ABSTRACT 


An electron paramagnetic resonance study of the *CH 2 CD 2 
radical in irradiated zinc acetate dihydpate has been made. The 
configuration of the radical in the crystalline lattice has been 
determined as nearly planar using a least squares fit. The HCH 
angle was found to be 116.9° + .3°. The principal elements of 
the oC-proton hyperfine tensors were calculated to be: 

-26.80 + .74, -59.59 + .44, -92.15 + .31 for -26.41 + .99, 
-58.98 + .45, -92.78 + .20 for H 2 . Three principal values of 
the g tensor were selected for the best least squares fit; 

2.00242, 2.00353, 2.0042. By experiment and an MO calculation 
the radical has been shown to be a IT-electron redical with the 
unpaired electron localized on the ot-carbon. Satellite lines 
were observed and interpreted as possible "spip-flip" transitions. 

A temperature dependence of the coupling constants for the 2 
protons was discovered from spectre taken at room temperature and 
at 133°K. Single crystals o^ Zn(CH2DC02)2‘2H20 were irradiated 
and an analysis of the spectra shoved that the barrier to formation 
of the radical vas the freaking of the C-’H bond on the undamaged 
acetate ion. 
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1, ESR Thi‘ory - Single Crystal Approach 
l.l Review of Radiation Damage Studies 

Researchers have found that Irradiation of certain 
materials creates one or more paramagnetic species whose study 
Is initially feasible by ESR methods (12), Experiments carried 
out with polycrystalline samples have unfortunately yielded 
rather limited Information on the hyperflne Interactions, l.e>, 
those interactions between nuclei and unpaired electrons. The 
broadening effect of the resonance lines due to the anisotropy 
present Is known to obscure the weaker signals (4), Single 

I i ‘ , 

crystals offer the best means of characterI*Ing these hyperflne 

. . . 
splittings In that radicals produced by Irradiation of single 

crystals have well-defined orientations In the crystal lattice. 

The external magnetic field can. then be located with respect to 

the radical. 

Since the first application of ESR Spectroscopy in 

I 

radiation damage studies about 14 years ago, work with single 
crystals have Included In part the saturated dlcarboxyllc acids 
and their salts, the simpler amino acids, and certain imsaturatod 
acids (6). The radical may be generated by Irradiation at room 
temperature or low temperature* and stability may be dependent 
on such factors as the envlrorutient In which the radical was created. 
Furthermore a radical generated by low temperature Irradiation doeii 
not necessarily mean that subsequent warm up will result in detecLion 
of the radical produced by room temperature Irradiation. The study 
of two radleaLs in irradiated crystals may be enhanced by the iact 
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that although both radicals are generated at room temperature, 
only one of the radicals Is present after low tmperature 
Irradiation. 

Once a radical Is generated In a known crystal lattice 

ESR techniques can extract Important Information about the radical. 

The measurement of the anisotropic hyperfine Interactions In the 

experimental spectra leads to the determination of one or more 

coupling tensors. The coupling tensor of the radical •CHCCO^H)^ 

In Irradiated malonlc acid was obtained from measurements of the 

splittings of resonance lines at various magnetic field orientations 

(11). The tensor's principal values,which are located along an 

orthogonal set of molecular axes In the radical, are related to 

the coupling tensor by a similarity transformation. Thus the 

orientation of the radical In the crystal lattice with respect 

to the molecular axis system Is partially known. The g tensor 

can be obtained by comparing the radical's hyperflne splittings 

with those of a standard whose ESR spectrum Is observed at the 
14 

same time. The N splittings of the three lines of peroxylamlne 
dlsulfonate In aqueous solution are used for such purposes (14). 

It should be mentioned that ESR studies of single crystals at 
controlled temperatures, after Irradiation at low temperatures, 
can yield Important Information about radiation damage processes. 

For further discussion of ESR methods In radiation damage studies 
Morton has an excellent review on the subject (12). 
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1.2 Problem Formulation using Spin Hamiltonian 

In dealing with a free radical in which a single 
unpaired electron may interact with several nuclei the effective 
spin Hamiltonian (using the notation of Slichter) is 

>1= (SH'yS - £ ^ S*Ac ‘Xi (21) 

L ‘ L 

This expression describes the interaction of nuclei and an 
unpaired electron in an external magnetic field for the case of 
quenched orbital angular momentum. The first term represents 
the Zeeman energy, the energy of interaction of th^ electron spin 
moment-With the external magnetic field. |3 is the Bohr magneton; 

S is the electronic spin; and g is the spectroscopic splitting 
tensor. .The second term is a summation of the Zeeman energies 
of the nuclei where is the spin of the ith nucleus and is 
the gyromagnetic.ratio of the ith nucleus. The third term 
represents the summation of the energy of the hyperfine interaction 
between the> electronic spin S and the spin of the ith nucleus 
is the hyperfine interaction tensor. The last term is better 
know as the combined effect of the isotropic Fermi contact inter¬ 
action and the anisotropic dipolar interaction. 

To obtain information about the A and g tensors experimental 
measurements,must be made on a single crystal in which the radical 
in this case has been formed by irradiation. Frpm the outward 
appearance of the crystal a set of orthogonal axes can be chosen 
(x, y, z) arid the crystal can then be rotated with knowledge of 
the position of the external magnetic field with respect to these 


11 








axes. It is possible to derive values of the elements of the A 
tensor in terms of the x, y, z system of axes (3), The trans¬ 
formation which diagonalizes this tensor will then define the 
relative orientation of the x, y, z system with jrespect to the 
radical principal orthogonal axes (X, Y, Z), Put one cannot 
determine from the ESR spectrum alone which of the principal 
axes is X, y, or Z. If the crystal struotpre of the molecule 
is known and it is assumed that the orientation of the atoms in 
the radical does not differ with that of the undamaged molecule, 
then these principal axes can be assigned. The investigation 
of the radical HCNH2"^C02 in irradiated glycine hy Ghosh and 
Whiffen was greatly enhanced by the known crystal structure of 
glycine (8). As in most radicals of this type the cayboxyl 
group in the undamaged glycine is assumed to be unmoved following 
irradiation (2). This is not an unlikely possibility and is a 
good basis on which to begin calculations of the direction cosines 
of the radical's principal axes with respect to the crystal axis 
system. These are the direction cosines that would be obtained 
from diagonalizing the A tensor. The principal values of the A 
tensor will lie along the principal axes in the radical. For the 
case of ano<.-proton coupling tensor the smallest value of the 
diagonal tensor corresponds to the direction of the C-H bond. 

Thus when the external magnetic field is placed along the C-H 
bond the interaction between the dipoles of the electron and 
nucleus is a minimum. If there are two nuclei with which the 
single unpaired electfon can interact then there will exist 2 
sets of principal axes, one for each coupling tensor. 
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The g tensor is anisotropic when the electron possesses both 
spin and orbital angular momentum (3). Thus in the electronic 
Zeeman interaction 

S does not represent the true spin of the electron since g tensors 
are anisotropic because the spin is not oriented exactly along the 
field direction. The true Hamiltonian for the magnetic interaction 
is 

BH-t + 

j I . 

where L is the orbital angular momenttan. The effects of this 
perturbing Hamiltonian may be represented by an "effective spin 
Hamiltonian" where S is considered as a fictitious spin. An accepted 
treatment is as follows: 

"W = Heff • S 

A 

Now the electron spin will be quantized along the unit vector h 

where Hg' is the magnitude of H and g' describes the effective 

magnetic moment of the molecule. Using the 3 principal values of 

the g tensor (g , g , g ) and the direction cosines of the 
XX yy zz 

T - 

magnetic field in the principal axis system (ot^, 

b, ^ 

The principal axes of the g tensor wilt not necessarily 
coincide with those of the hyperfine tensor. 

Once the g and A tensors are known the effective spin 
Hamiltonian can be used to calculate the probable transitions 
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that will occur and thus the spectrum for a particular orientation 
of the external magnetic field can be predicted* 

1.3 Solution of the Spin Hamiltonian 

With the A and g tensors for a particular molecule, 
a method must be used to obtain the transition frequencies and 
their probabilities between various energy levels in the molecule. 

r 

One such method can be outlined as follows: 

a. Consider the Interaction of a single 
unpaired electron with n nuclei and form 
the basis functions to be used In setting 

up the Hermltlan matrix. The basis functions 
will be of the form 

. 

l.e., the product of the electrpn ppln function 
and n nuclear spin functions. The matrix will 
be a square matrix with 2 elements on a side 
where x = nvtnber of spin functions. 

b. Operate on the terms in the effective spin 
Hamlltpnlan with all required permutations 

of the basis functions. The diagonal elements 
are 

. 

i. 
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The ofI diagonal elemonts will contain 
Imaginary elements due to the Hermltian 
properties of the system. The Hermitlan 
matrix can be diagonalized by a method 
described by Wilkinson (33). Using this 
method a real matrix is created and the 
resulting eigenvalues and eigenvectors 
occur in pairs. 

c. Set up the transition probability matrix 
and transpose to the diagonal basis. 

Since the eigenvectors are Imaginary 
the transformation can be handled as 
follows: 

Eigenvector » U +• LV 

Transition probability matrix ■ P 

d. Square the elements of the transition 

probability matrix in the diagonal basis. 
The absolute value of squared is the 
relative intensity of the transition 
occurring at frequency A] •* • 

The Alsare the eigenvalues expressed in 
megacycles. With frequencies and 
intensities the spectral envelope can be 
drawn using the derivative of a Lorentalan 
distribution. 
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The following alternative solution of the Spin Hamiltonian 


was used by Box (2) In his treatment of deuterated glycine, 


'*ND^CH 2 C 00 ". The effective spin Hamiltonian Is the same as described 


In section 1.2. By Ignoring the components of S perpendicular to 
H the energy levels of the spin system are given by 


E = P H % likl ktl 


where g Is the observed g value, Is the component of S along H, 


and Is the component of along h^ where h^ Is defined by the 
equation 



A, B and C are the principal values of the hyperflne tensor, A^. 
The unit vectors u^, v^, w^ are the vectors along which each of 
the principal values A, B, and C lie In the crystal axis system. 
The direction cosines of the external field with respect to the 
crystal axes are given by o^^, |3^, The relative probabilities 

of the transitions are functions of the angle between the vectors 
h^ (+) and h^ (-), the values of h^ corresponding to equal to 
+ \ and - %, respectively. 

A computer program was developed using the above theory 
to calculate the relative Intensities and frequencies for the 
transitions which occur (23). Interactions are considered between 
a single unpaired electron and any number of nuclei. However, the 
orientation of the radical In the crystal must be known In order 
to calculate the vectors (+) and h^^ (-). The unit vectors 
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^1’ * '*^1 direction coginea which pli^ce the principal 

values of the A tensor along the crystal axes. The program, 
shown in Appendix 1, does not calculate the g anisotropic effect 
but assumes g « 2.0023, Hereafter, this program shall be referred 
to as moiCT. 
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2. Background on Zinc Acetate Dihydrate 

Many investigations have been focused on the hyperflne 
interaction between o(. - protons and a single unpaired electron 
in the 2 p - orbital of the central carbon atom for radicals 
such as HC-CCOOH)^ (11) • In many of the systems studied several 
IT -electron radicals were created by irradiation thus making 
Identification of spectra difficult. The ‘CH^CO^ radical has 
been observed with the NH^'^CH^CO^ radical in irradiated glycine 
(13). In working with the glycine single crystals, Morton 
concluded that the two oC-proton coupling tensors of .CH 2 COO” 
have axial symmetiry and suggested that the CH 2 group was 
experiencing hindered rotation about the C-C bond. Box and 
Freund demonstrated an application of ESR theory to the radical 
in irradiated glycine and had success in assuming the radical 
was planar. Horsfield and co-workers showed that for the 
•CH 2 CO 2 H radical in maIonic acid, the plane of the CH 2 group 
was approximately perpendicular to the plane of the (OO 2 H) group. 
A value of 116.5^+5° was obtained for the HCH angle in the 
‘CH2C02H radical (9). The principal values of the 2 coupling 
tensors obtained for this radical were 


^1 

»2 

•30 

-37 

•55 

-59 

■91 

CM 

a\ 

1 


Morton calculated for the CH 2 CO 2 radical (planar) the following 
principal values for the coupling tensors, assuming the tensors 
to be equivalent (13). 
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and 


-48 

-56 

t 

-69 

But Box and Freund successfully predicted spectra for the 
planar •CH^CO^ in glycine by'using the principal values of the 
coupling tensor obtained by McConnell (11) for *011(00011)2. 

-29 

1 

-61 

-91 

Valenti discovered that the CH^CO” radical existed in 
irradiated Zinc Acetate Dihydrate (28). The radical was found 
to be stable at room and liquid nitrogen temperatures. The 
spectra taken as the temperature wa? lowered showed the center 
line of the 1:2:1 hyperfine pattern splitting into 2 lines. 

This confirmed previous observations that a barrier to internal 
rotation must exist (13), The use of 1st order perturbation 
theory yielded the principal values of the two c<-prQton 
coupling tensors but the spectra was tpo complex for Q^ny further 

.9 . 

analysis with this theory. This was the motivation for a complete 
investigation of the CH 2 CO 2 radical, but specifically: (a) to 
explain the suspected anistropic effects in the complex apectra 
for various orientations of the external magnetic field; (b) to 
predict by a theoretical calculation the exact spp^ctra of the 
CHzCO” radical; (c) to obtain not only the best principal values 
of the A and g tensors but the correct orientation of the radical 







in the crystal within the error limits of a least squares fit; 
and (d) then to prepare ZnCCH^I^CO^)^crystals for exchange 
rate and barrier studies. 

The orientation of the radical in the crystal was calculated 
by starting with the information given by van Niekirk, et. al (18). 
The crystal structure of ZnCCH^CO^)^is monoclinic with space 
group C 2/c and unit cell dimensions of a = 14.5, b - 5.32, c = 11.02 
angstroms with P = 100° (Figure 1). The Zn atoms lie on a 2-fold 
axis and although there are 4 formula units per unit cell, only 
two magnetically inequivalent sites exist. The relationship of 
the orthogonal axis system to the cell axis system is 

a = sin^ X + cos 3 fc 
b = y 


X will also be referred to as a*. With these relationships the 
coordinates of the atoms in the x, y, z system were calculated 
to be 



X 

y 

z 

Zn 

0.0 

.665 

2.755 


1.014 

1.064 

.868 


1.471 

2.250 

2.683 


1.613 

2.128 

1.314 

Co 

2.614 

2.804 

.454 
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A 

b 



Figure 1, The orthogonal crystal axis system (lyz or a*bc) 
and Its relationship to the crystalline axis 
system (ahc)* 
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The vectors fpr the distances between these atoms were 
calculated and when normalized became direction cosines of 
the vector with respect to the crystal axis system. Thus for 
example: 

Vc^c^ = (.6751, .4559, -.5800) 

Vo^02 = (.2062, .5352, .8191) 

The angle between the 

(Figure 2). At this point the assumption was made that Vc^C 2 
is unchanged during radiation damage to the crystal. Then the 
HCH plane was arbitrarily set up to include the 
and to be parallel to Vo^- 02 . The planar configuration was 
chosen since Valenti could not explain his results using the 
perpendicular model (28). The perpendicular model was not 
ruled out. 

A vector perpendicular to the CH^ plane is 
= V0j^02 X Vc^c^ = (-.6878, .6755, -.2685) 

The directions of the C-H bonds were established assuming 
an HCH angle of 120°. 

Vp = Vcj^c^ X Vj, = (.2694, .5802, .7696) 

Vc^H^ = ^ (-.1043, .2745, .9565) 

VC 1 H 2 = -*5^0^02 - 3/2 Vp = (-.5709, -.7305, -.3765) 

Since the principal values of an oC-proton cpupling tensor change 
very little from one saturated hydrocarbon to another (if delocalis¬ 
ation of the odd electron is not critical), those values measured 
in the malonic acid radical, •CH(C 00 H) 2 > were used. 
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Figure 2. Locations of the hyperflne axes in the *CHoC0 
radical(planar configuration), ^ 
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A = -29. 


B = -61. 

C = -91. 

Figure 2 shows the CH 2 CO 2 radical and the assumed orientation 
of the molecule with the principal axes assigned as shown. 

The principal values A, B, C, were designated to be along the 
C-H bond, the p^ orbital axis, and in the CH 2 plane, perpendicular 
to the C-H bonds. These directions correspond to vectors u, v, 
and w in the crystal axis system. Thus the direction cosines 
for the A tensor were calculated to be 




X 

y 

z 

«1 


-.1043 

.2745 

.9565 


^1 

-.6878 

.6755 

-,2685 



.7194 

.6849 

-.1175 

«2 

"2 

-.5709 

-.7305 

-.3765 


^2 

-.6878 

.6755 

-.2685 


”2 

.4500 

.1047 

00 

00 

• 

t 
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3, Experimental 

3.1 Growing of crystals 

Various methods were used to grow the crystals needed 
for use in the spectrometer. All of the methods depended on 
the slow evaporation of water from saturated aqueous solutions. 

A 2.5:1 mass ratio of the zinc acetate to water worked best. 

A 115 ml solution in a 600 ml beaker resulted in well formed 
crystals but required from 2-7 days. If watch glasses are 
used, crystals can be gr'owu overnight easily. Usually the 
above procedures are successful only on the first cycle because 
pf the contamination of the solution on standing. 

’ Another proven method was adding ZnO to 3-3.5 g con¬ 
centrated acetic acid in 80-100 ml until a neutral solution 
resulted (no excess ZnO should be added). The solution was 
allowed to stand in a 150 ml beaker for 48-72 houps with the 
”hood** sliding window 1 inch above the beaker. Perfect crystals 
were produced evefy time this method was tried. 

Acetic anhydride and ZnO were used together with D^O 
to obtain the crystals for study of the effect on the satellite 
lines observed in the spectra. The crystals (Zn (CHoCO^) *20 0) 
were grown by slow evaporation of the solution frffn watch glasses. 

3.2 Identification of Crystal Axes 
3.2.1 X-Ray Diffraction 

The putward appearance of the monoclinic crystals 
of zinc acetate dihydrate showed clearly which axis was the b axis, 
the two-fold axis. But to insure that this identification was 
correct an x-ray diffraction pattern was taken on an oscillation 


25 









photograph. The crystal was glued to a glass fiber so that 
the incident x-ray beam was perpendicular fo the assumed b axis. 

The crystal was rotated about the b axis and a photograph was 
taken. From the distances between the layer lines in the photo¬ 
graph the repeat distance of the Bravais lattice along the 
direction of the b axis was obtained. By a simple calculation 
using Braggs Law the length of the unit cell along the b axis 
was found to be 5.18A°, which compared favorably with the value 
(5.32A^) found by van Niekirk, et. al., for the unit cell dimension 
along the b axis (18). 

3.2.2 Crystal Alignment 

The orientation of the crystal in the spectro¬ 
meter cavity was critical in order to determine the precise 
location of the magnetic field with respect to the crystal axis 
system. Two angles of rotation were needed to completely describe 
the ^-dimensional system. The crystal was rotated about one angle 
in the cavity but the other angle was varied outside the cavity. 

A Baush and Lomb microscope was set up with a revolving platform 
as a work pedestal. The turntable was calibrated such that 1 
revolution = 360 units or degrees. The crystal was mounted on a 
quartz rod (28) with the (100) face parallel to the axis of the 
rod. The rod was clamped to the pedestal so that the crystal was 
centered on the turntable and the rod itself was positioned 
perpendicular to a tangent on the periphery of the turntable. 

The crystal was oriented on the rod by locating the b axis vith 
respect to the axis of the rod with the aid of an eyepiece with 
crosshairs. The axis of the rod was positioned parallel to one 


26 




of the crosshairs♦ The turntable was rotated the desired amount 
in degrees and the edges of the crystal parallel t^o the b axis 
were lined up with the reference crosshair. 

3.3 Radiation Dosage 

The single crystals were exposed to x-irradiation at 
50KV and 20 ma from 2-4 hours using a molybdenum target. The 
crystals were placed on scotch tape and mounted in the open windows 
of the x-ray machine allowing the port covers to be closed at all 
times. The irradiations were made at room temperature. 

3.4 ESR Spectrometer 

The ESR spectrometer used was a Varian V-4502-13 with 
aFleldial V-Fr 2503 magnetic field regulation and a nine inch 
magnet, using 100 KC modulation, and a conventional reflection 
type cavity equipped with a Varian variable temperature control 
unit. 

The quartz rods on which the crystals were mounted had 
two flats ground parallel to each other on both ends of the rod. 

The end opposite the mounted crystal was clipped to an angle 
indicator positioned over the center of the spectrometer cavity. 

It was estimated that the angular position of the crystal could 
be measured accurately within + 1° with this apparatus (28). 

The crystal was positioned in the yz plane by use of 
the microscope. The angle of rotation was designated as (ji) 
and the angle of rotation corresponding to the angle indicator 
above was called (&) (Figure 3). All spectra taken were labeled 
hy (^) for easy identification. The crystal was mounted with 
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Figure 3* Single crystal of Zn(CH3C02)2*2H20 mounted on a 
quartz rod with the h axis coincident with the 
axis of the rod. Crystal rotated by ^ and rod 
rotated by & • 
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its X axis extending into the rod and perpendicular to the axis 
of the rod. A complete data run was made by starting with the 
magnetic field along the -x axis((^) direction) and incrementing 
by 15^ until the field was along the -he axis - this was done for 
^= 0, 15, 30, 45, 60, 75, 90 degrees. Each spectrum was displayed 
on an x-y plotter and the magnetic field reading for each line 
was taken from a gauss field dial on the spectrometer. 

3.5 Synthesis of Zn (CH2DC02)2*2H20 

Ethylene oxide^CH^CH^ was chosen as one of the starting 
materials because its ring is easily opened by a suitable reducing 
agent (7). LiAlD^ was chosen as the reducing agent and the follow¬ 
ing reaction was carried‘out in a 500 ml, 3-necked flask immersed 
in an ice bath; 


4CH.-CH^ 

V 


+ L\h\0^ 


4CH2DCHgOH -I- sul-fft+es 


The solvent used was triglyme |J., 2-Bis (2-methoxyethoxy) ethane^ 
since its boiling point (220-225^0) was well above that of water 
and CH 2 DGH 2 OH. The solvent was purified by first using KOH pellets 
and letting solution stand for 72 hours. Approximately 250 ml 
was poured into a distillation flask and 6 g LiAlH^ was added. 

After the initial reaction of the H 2 O and LiAlH^ had subsided, 
the solution was distilled under vacuum to obtain the high boiling 
triglyme. The operating pressure was .43 mm Hg, The distilling 
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temperature was 61“62°C. 50 ml of the solvent was then placed 

in the 3-necked flask, with 3 g of LiAlD^, The ethylene oxide 
was added to 50 ml of chilled triglyme in a 250 ml round bottomed 
flask on a "top loader" balance so that the ethylene oxide could 
be added by weight. 12.5 g of ethylene oxide was added to the 
triglyme and the solution was poured into a dropping funnel on 
the reaction vessel. The solution was allowed to run into the 
flask slowly, A condenser was placed on the reaction vessel with 
a 0 °C water flow in order to minimize the loss of ethylene oxide 
(b.p.lO°C). After the reaction appeared complete, 58 ml of 
3M ^ 280 ^ was added and the Increase in viscosity that resulted 
indicated that the sulfates were coming out of solution. The 
CH 2 DCH 2 OH and water were distilled over at a temperature range 
of 78 to 102°. Yields were calculated to be around 96% giving 
approximately 13 g of CH 2 DCH 2 OH in the water alcohol mixture. 

The deuterated alcohol was then added slowly by burette 
to an oxidizing mixture of 

40.7 ml H 2 S 0 ^ 

56.2 g K2^r20^ 

212.0 ml H 2 O 

in a 3-necked flask immersed in an ice bath. The following 
reaction took place while keeping the solution in the vicinity 
of 25°C. 

^CH^DCHgOH + -h8H^S0^ = 

BCHjDCOOH + E K2 SO^ H HgO (iO) 
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The solution was refluxed for 5-10 minutes and then the CH 2 DCOOH 
and water were distilled over at 100°C. Yields ranged from 65 to 
75 % giving 10-12 g of the acid. 

ZnO was added to 85 ml of the acid solution (3-3.5gCH2DCOOH) 
until the solution turned neutral but not cloudy (excess ZnO). 

ZnO + ZCH^DCOOH = Zn CCH^DCO^.'^g^-I-HpO 

The slow evaporation of the product solution from a 150 ml 
beaker produced the needed crystals for ESR study. 

Vapor-phase chromatography was used to identify the alcohol 
and acetic acid in their product solutions and the peaks in comparison 
to those of a standard were indicative of relative amounts present 
also. The VPC analysis did not contain any impurity peaks. 
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FLAP 


Figure 4. Bending vibrations of the CH 2 group in the •OH 2 OO 2 
radical. 
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Problem Solution and Spectral Analysis 
4.1 Characterization of the '0112002 radical 
4.1.1 Theoretical Formulation 

The main computer program developed in this 
investigation was constructed around two main subprograms; 

(a) a method for prediction of the transition frequencies and 
intensities (PREDICT); and (b) a least squares method for 
obtaining the best fit between calculated and observed frequencies. 
This maiu program shall be hereafter referred to as OPTIM. The 
method for prediction of transition frequencies and intensities 
was available in PREDICT and the input to this program required 
a precise orientation of the CH 2 group, A change in this orien- 
tation produced a change in the predicted transition frequencies 
and intensities because of a simultaneous change in the direction 
cosines. Thus a means of describing the location of the CH 2 group 
was needed in order to perform any repetitive calculation such as 
a least squares fit. The use of a least squares fit required that 
certain parameters be chosen and the angles concerned with the 
location of the CH 2 group were considered good choices. The selection 
of parameters was based on the criteria of at least determining 
whether the CH 2 group was planar or perpendicular to the rest of 
the molecule. The twist of the CH 2 group was then chosen as the 
first parameter which involved a rotation about the C-C bond. 

(Figure 4) Three more parameters were chosen to account for the 
wag, flap, and HCH angle of the CH 2 group. The relationship of 
these three parameters to the direction cosines of the CH 2 principal 
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axes 


were solved by using a similarity transformation technique. 
Any random vector (a, b, c) can be rotated to a position along 
the z axis by the following transformation 



The following matrix multiplications illustrate the 
a transformation matrix which will rotate any other 
number of degrees 9 about the vector (a, b, c). 





be 

a 



formation of 
vector a 


d 

hSL 

a 

b 


0 


\ 


-d 



This transformation was permitted to operate on the direction 
cosines of the g and A tensors for each value of twist, wag, and 
flap. The direction cosines were updated each time a new set of 
parameters was calculated. The HCH angle was not involved in a 


rotation but was used to set_ up the planar configuration for 

each new set of parameters. This configuration was assumed because 

initial calculations using PREDICT showed the planar form to yield 

the more correct "stick spectrum". 

PREDICT was modified to include the calculation 

of an effective g value which was necessary due to the 2 magnetically 

inequivalent sites of Zn(CH 3 C 02 ) 2 * 2 H 20 for orientations of the 
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magnetic field in other than the xz plane. The transitional 
frequencies computed by PREDICT for each site were shifted 
upfield or downfield depending on the effective g value. The 
calculation of an effective g value was made in section 1.2. 

The principal values of the g tensor used as initial parameters 
in this investigation were 2.0020, 2.0042, and 2.0034, which 
were values corresponding to directions perpendicular to the 
CH^ group, parallel to the C-C bond, and in the plane of the 
CH 2 group perpendicular to the C-C bond. Due to their import¬ 
ance these variables were also selected as possible parameters 
along with the 6 principal values of the 2 coupling tensors. 

Both coupling tensors were included because other work has shown 
that hyperfine coupling with -protons is highly anisotropic (16), 
but the principal values are expected to be similar (3). 

The transitional frequencies computed by PREDICT 
were placed in a matrix to be used to calculate differences between 
selected frequencies. These differences were then compared with 
observed differences taken from the experimental spectra and the 
results formed a reference matrix to be used in the least squares 
fit, which was developed for use with linear or non linear equations. 
In this investigation a highly non .linear situation was assumed to 
exist and the derivative matrix was formed by varying each parameter 
only a small amount. The mathematical formulation of the least 
squares method is as follows; 
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i 


y 

y 


i 

o 

i 



= error between calculated and observed value 
= ith calculated point 
= 1th point, observed value 
= jth independent parameter 


= . / c)a . 

1 J 



-ith point 



The least squares restriction: 







+.idiidcKAaj=0 

Define D = matrix of elements 

E = column matrix of elements 

= column matrix of elements ^a^ 

D*E + D*D^x « 0 
Thus/\x = - (D'D)‘^n'E 

The matrix Ax represents the best change in each independent parameter 
for the best fit of the data. 

The elements of this matrix were added to the corresponding original 
parameters and the entire program of OPTIM was repeated until successive 
iterations produced only very small changes in the parameters. The 
derivative matrix (d ) was generated by introducing incremental 

1 J 

changes in each parameter while holding others constant and calculating 
the resulting effect on the differences between the transition frequencies. 
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Wulff Net 



Figure 5. Stereographic net used to catalogue magnetic 
field orientations. 
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A measure of the uncertainty In each parameter 
can be calculated using a probabalistic approach. From the 
least squares analysis the best changes In each parameter for 
the best fit were 

Ax = (D'D)‘’'D'E 

If there were no error Involved 

where E * calc -obs points 
but, in general, E ■ D^x + e 

Uncertainty in Ax s<.Ax'> 

"True" value of AX — AX 

A 

Value of Ax at some time is AkX 

- Average of (Ak - (^Ax - 

- Average of jCD'D)"^ D'EI -Ax] D'E -Ai^' 

“ (D^D) ^ ©"P ea')D(D*Dy^ 

- cr^ (D'D)'^ 

2 

where<r is the standard deviation squared (avg of ee'). This 
expression shows that if D'D is diagonal then the errors are truly 
determined. 

4.1.2 Data Obtained 

With the crystal mounted as described In section 
3.4, ESR spectra were taken at 15° Increments of @ and 0. For a 
corresponding value of 9, 0 was varied from 0 ° to 180° and the 
resulting number of possible magnetic field orientations was 79. 

The number of spectra was so extensive that a stereographie net was 
used to catalogue the different orientations ( 0 ). A typical Wulft 
net Is shown in Figure 5 with the crystal axis system labeled and 
some of the field orientations indicated. 


38 





TABLE I. Data used for the least squares fit 


Field 

Orientation 

es) 






Field Dial 

Differences 

Field 

Field Dial 

Differer 

Reading* 

Taken 

Orientation 

Reading* 

Taken 

20.50s 

29.78 

(T) 

17.84s 

20.96 

30.84s 

10.34 

38.80s 

26.61 

50.28s 



44.45s 


60.61s 



65.38s 


22.55s 

23.50 

(S 

14.19s 

18,50 

27.20W 

17.17 

19.70W 

28.62 

35.10W 


37.24s 


46.05s 



38.30w 


50.60W 



43.40s 


58.50W 



48.30W 


63.20s 



56.95W 

61.84s 


22.20s 

27.45W 

40.58 

/15\ 

13.50s 

19.39 

36.80W 


\90j 

17.40s 

15.49 

42.49s 



27,80w 

29.21 

48.12W 



32.89s 


57.20W 



37.88s(2) 


62.78s 



42.81s 

62.00s(2) 


21.60s 

17.70 




26.01W 

24.09 

(30\ 

12.54s 

18.31 

34.29W 


\90j 

17.98s 

31.75 

39.30s 



25.78m 

12.87 

45.69s 



30.85s 

26.31 

50.23W 



39.49s(2) 


48.68w 



44.29s 


63.29s 



49.26W 

62.10s(2) 


20.16s 

18.12 




24.95W 

26.07 

(U5\ 

12.08s 

17.58 

33.19W 


[ 90 ) 

18.43s 

32.69 

38.28s 



24.84m 

11.23 

46.23s 



29.66s 

26.34 

51.02W 



34.71m 


59.38W 



39.61s(2) 


64.22s 



44.77s 

49.54m 

61.83s(2) 



*Gauss field dial on spectrometer. 

Readings taken from a 100 gauss sweep. 

s, strong; w, weak; ( ), no. of lines overlapping. 
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Field 

Orientation 

Field Dial 
Reading* 

Differences 

Taken 

Field 

Orientation 

Field Dial 
Reading* 

Differences 

Taken 

/60\ 

11.90s 

17.42 

/120\ 

30.85s(2) 

10.02 

V90/ 

17.40s 

32.71 

Vl5 J 

46.40s 



28.32s 

11.92 

y 

51.25s(2) 



38.90s(2) 

27.21 


56.42s 



44.61s 



72.00s(2) 








61.25s(2) 

12.40s 

16.09s 

24.46W 

30.02s 

35.46s 

41.09s 

56.03W 

61.00s(2) 

22.80s 

26.95s 

46.20s(2) 

54.60s(2) 

69.30W 

75.65s(2) 

23.20s 

28.00s 

43.68s(2) 

50.63W 

56.00s(2) 

72.45s 

76.20s 

25.10s 

30.11s 

42.25s(2) 

57.99s)[2) 

70.53s 

75.08s 

30.11s(2) 
43.27s 
47.04w 
58.47s 
71.73s(2) 


17.80 

19.37 

14.11 

(") 

26.19W 

31.39s(2) 

47.08s 

51,40s 

56.30s 

68.66w 
73.00s(2) 

41.61 

4.15 

(IS 

25.51W 

30.20s(2) 

48.63s(2) 

56.30s(2) 

69.00W 

73.95s(2) 

43.75 

53.00 

44.45 

48.20 

(s) 

21.70s 

28.18s 

47.03s(2) 

51.35s 

55.25s 

68.60w 

73.60s 

77.63s 

55.93 

51.90 

49.45 

49.98 

40.42 

45.43 

44.97 

(s) 

21.58s 

29.63s 

45.26s(2) 

53.40s 

56.78s 

72.258 

77.148 

55.56 

42.62 

50.67 

47.51 

15.20 

&) 

23.37s 

30.56s 

44.63s(2) 
55.78s(2) 
71.40s 
76.75s 

53.38 

40.84 

48.03 
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Field 

Orientation 

Field Dial 
Reading* 

Differences 

Taken 

Field 

Orientation 

Field Dial 
Reading* 

Differences 

Taken 

fa 

16.30s 

12.15 

/130\ 

21.00s(2) 

25.50 

19.80s 


\ 45 7 

34,60s 

26.00 

33.25s 



39.40s 

30.30 


39.25s(2) 



46.48s(2) 



45.40s 



60.60s 



58.80s 



64.90s 



62.30s 


fl35\ 






14.80W 

16.35 

(150\ 

12.90W 

12.08 


21.78s(2) 

30.30 

V30/ 

18.60s(2) 

0.00 

33.45s 



33.32s 



40.97s 



39.52s(2) 



44.90s' 



45.40s 



49.80s 



60.40s(2) 



60.13s 






63.75s 


f45\ 

21.40s 

56.79 

ri50\ 




29.80s 

25.05 

13.13w 

14.57 


49.25s(3) 

43.80 

\45/ 

19.50s(2) 



54.85s 

33.45 

32.86s 



73.60s 

52.20 


41.19s(2) 



78.19s 

43.39 


47.43s 




23.34 


61.87s(2) 






67.30W 


feo] 

21.32s 

57.28 





30.95s 

24.80 

fl5\ 

26.05s 

21.50 


47.84s(3) 

41,95 

UoJ 

29.36s 



55.75s 

34.41 


48.05s 



72.90s 

9.63 


53.10s 



78.60s 

51.58 


76.30s(2) 




47.65 






22.87 

/30\ 

'23.50s 

23.27 




Uoj 

30.51s 

16.91 


22.76s 

55,54 


47.42s 

7,01 

UsJ 

31.60s 

41,10 


53.58s(3) 



48.30s 

46.90 


74.93s 



53.92s(2) 

8.64 


78.00 



72.70s 

49.74 

f45\ 




78.50s 


22.10s 

25.68 




Uo/ 

31.37s 

56.45 

17.90s 

48.27 


47.78s 

42.58 

32.90W 

27.70 


52.82s(3) 

47.18 


38.38s(2) 



73.95s 

51.85 


45.60s(2) 



78.55s 

9,27 


61.10s 

66.17s 
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Field 

Orientation 

~S“ 








Field Dial Differences 
Reading* Taken 


21.60^- 

25.85 

31.42^- 

56.35 

48.54s 

41.85 

52.05s(3) 

46.98 

73.08J 

s 

51.48 

78.40s 

& 

9.82 

22.15s 

5 

55.10 

31.10s 


41.32 

49.85s(3) 

46.56 

72.08s 

3 

49.63 

77.56s 

B 

8.85 


17.40s 

49.08 

21.85s 

11.80 

33.65s 

32.83 

41.80s(2) 


49.30s 


61.00s 


66.48s 



20.90s(2) 

19.38 

31.62s 

29.48 

42.05s(2) 

34.04 

51.00s 

14.66 

61.108 


65.668 


19.858(2) 

19.70 

30.10s 


41.478(2) 


49.80s 


62.358 


68.50W 



Field Field Dial 

Orientation Reading* 

TioV 23.55s 

1,75 j 30.30 b 

45.00s 
51.908(2) 
74.59s(2) 

MSX 21.95s 

^75/ 30.70s 

43.90s 
51.82s(2) 
72.72s 
76.82s 


21.408 
\75/ 30.28s 

44.218 
52.048(3) 
67.45W 
72.188 
76.658 


( 75 \ 23.278 

V75/ 30.968 

47.488 
54.408(2) 
59.SOW 
72.528 
. 76.75s 

/ 120 \ 19 . 008 ( 2 ) 

\ 75 ) 29.75s 

41.388(2) 

49.508 

63.90i 


18.328 15.63 

21.30s 

32.60& 

42.30s(2) 

48.23s 

63.l5s(2) 

25.80s 21.06 

29.80s 16.92 

41.10W 4.08 

46.728 
52.26s(2) 

55.23s 
75.80s(2) 


/135^ lO.lOw 

[ 75 J 18.90s 

29.138 
40.678(2) 
47.55s 
61.60s(2) 

/l50\ ll.lOw 

\ 75 / 16.288 

20,408 

31.808 

41.47s(2) 

47.308 

62.758(2) 

67.958 
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Differences 

Taken 

21.83 

14.70 

6.75 


21.95 

54.87 

42.02 

13.20 

8.75 

50.77 

46.12 


22.81 

55.25 

41.90 

13.93 

8.88 

50.78 

46.37 


53.48 

42.16 

46.39 

49.25 


19.75 

34.15 


18.42 


15.52 

31.02 

11.40 









Field Field Dial Differences Field Field Dial 

Orientation Reading* Taken Orientation Reading* 







27.48s 

46.04 

31.30s 

38.65 

42.20s(2) 

42.22 

58.75s(2) 

42.47- 

69.95s 


73.52s 


26.30s 

49:90 

31.20s 

40.40 

45.00s(2) 

45.00 

57.07s(2) 

45.30 

71.60s 


76.20s 


25.20s 

52,65 

31.86s 

40.64 

47.60s(2) 

45.99 

53.90s(2) 

47.30 

72.50s 


77.85s 


24.42s 

52.48 

30.88s 

40.96 

42.03w 

46.02 

46.52w 

47.42 

50.52s 


55.67w 


71.84s 


76.90s 


25.67s 

51.16 

30.22s 

42.59 

45.55s 

46.61 

48.85s 

47.14 

53.87s 


57.77s 


72.81s 


76.83s 
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Differences 

Taken 


in 











Each spectrum was analyzed with the aid of the 
"stick spectrum" predicted by PREDICT for a particular field 
orientation. The "stick spectrum" was plotted after arbitrarily 
assigning sites 1 and 2. The best lines were selected for 
taking differences; some spectra were not used because of too 
much overlap. Differences were taken within each site and between 
the 2 sites. These differences are tabulated in Table I for each 
spectrum that was usable in obtaining the least squares fit. The 
gauss field dial readings (from spectrometer) taken for each 
experimental spectrum are also tabulated. Th^se readings correspond 
to the relative positions of resolvable lines in a 100 gauss sweep. 
The experimental spectra were taken at 133°K at a magnetic field 
of 3210 + 50 gauss. 

4.1.3 Correlation of data with Theory 

The experimental spectra for selected field 
orientations are shown in Figures 7-9. The predicted spectra 
are plotted also to show the close fit that was obtained by 
least squaring the differences taken. Although 3 iterations 
were used, the most significant convergence of OPTIM occurred 
in the first iteration. The best fit was obtained using 166 
differences selected from 55 different field orientations. 

More than one crystal is represented by these data and it was 
significant that two different irradiated crystals did not 
differ in their spectra and the differences obtained were with¬ 
in + .5 gauss of each other. After 3 iterations of the least 
squares program the matrix containing the calculated minus 
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Wulff Net 



Figure 


Errors 

protons 


In. coupling constants for the meth 7 leno 
^gjln *012002 for a 1° error In angle p . 
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observed differences 

was as 

follows: 

the error differences 


are read by 

rows and 

correspond to the same 

order as the 


observed differences 

In Table I. 





.35 

-.01 

.26 

-.48 

.25 

.12 

-.21 

-.10 

-.09 

-.28 

.22 

.30 

.77 

-.49 

.28 

.13 

-.19 

-.34 

.11 

.29 

-.22 

-.03 

-.15 

.08 

-.22 

.01 

-.32 

-.06 

.26 

-.37 

-.45 

.30 

-.30 

-.56 

.57 

.34 

-.33 

.23 

-.35 

-.32 

.04 

.27 

.03 

.39 

-.09 

.22 

-.04 

-.08 

.26 

.14 

-.30 

.04 

.12 

.28 

.84 

.31 

.14 

-.09 

.45 

.32 

.22 

-.00 

-.14 

.09 

-.55 

-.20 

-.44 

.09 

-.11 

-.00 

.53 

.20 

-.37 

-.03 

.23 

-.14 

.25 

.22 

-.04 

.61 

-.02 

-.11 

.53 

.30 

.49 

.38 

-.14 

-.13 

-.47 

.00 

.11 

.02 

.10 

-.01 

.48 

.56 

-.26 

.07 

-.03 

.04 

.53 

-.13 

.35 

-.03 

-.13 

-.04 

.87 

-.14 

.28 

.26 

- ,06 

-.34 

-.50 

.26 

-.21 

.34 

-.49 

-.64 

-.26 

.00 

-.33 

-.64 

.23 

-.10 

-.23 

.00 

-.41 

-.78 

-.88 

.16 

-.33 

-.45 

-.29 

-.43 

-.31 

.38 

.02 

.05 

-.68 

-.44 

-.11 

.71 

.46 

.24 

-.48 

.82 

.04 

.31 

.07 

-.48 

-.50 

.10 

-.10 

-.36 

-.09 

-.36 

.34 

-.08 

.38 

-.12 

-.36 

.50 

.07 

.06 





The standard 

deviation for 

this data 

was calculated to 

be .335 

gaus s. 

The angle 0 Involved 

In the 

rotation 

of the 

crystal In 

the spectro- 


meter and defining the center of a transition were considered the main 


sources of error. An error analysis using the proton coupling constants 
Is shown In Figure 6. The errors In the coupling constants for protons 
1 and 2 for a 1° error In angle 0 are displayed on a Wulff net for 
selected orientations. In general, as the angle varies, the error 
In one coupling constant Increases while the error In the other 
coupling constant decreases. An average error estimate for each 
coupling constant Is + .15 gauss. The coupling constants for field 
orientations close to the b and c axes are very sensitive to changes 
In 0, But as the field vector approaches either the -a'^ or a*^ 
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^ 7, ESR spccLra oi lIu’ radical at L33 K, Magnetic ILeUl perpendicular to two-Cold axis. 
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ESR .spectra of the -CIl CO' ratiical at 133®K. Magnetic field pcrptMidicular to z axis for /45\tfpectrnm 
and perpendicular to x axis Cpr other spectra. V90/ 
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Fifjtjrr* 9. ESR spectra of the •ClI^CO* radical at 133^K for selected magnetic field orientations. 
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direction the coupling constants are very insensitive to change. 
Here the field vector is entering the plane perpendicular to the 
b axis where both sites are equivalent. Spectra taken with the 
field close to or along the b axis were the most sensitive to an 
error in 0, The effect of an error in 0 on the coupling constants 
was found to be negligible. 

The error in defining the center of each trans¬ 
ition was estimated as + .20 gauss which makes the total error 
estimate = + .35 gauss. This is comparable to the standard 
deviation obtained in OPTIM. 


4.1.4 Results 

The parameters obtained using the data shown 
in Table I are listed below with the uncertainty calculated for 
each value: 

Orientation of the CH2C0" radical Uncertainty 


Twist 

- 5.25° 

.38 

Wag 

- 5.00° 

.18 

Flag 

- 9.88° 

.26 

/hch 

116.90° 

CM 

• 


Principal Values of; Tensor 


-26.80 

.74 

-59.59 

.44 

-92.15 

.31 

Tensor 


-26.41 

.99 

-58.98 

.45 

-92.78 

.20 
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Direction Cosines for; 


( 

.547, 

.722, 

.423) 

( 

.768, 

-.634, 

.088) 

(- 

.332, 

-.276, 

.902) 

( 

.048, 

-.081, 

-.996) 

( 

.768, 

-.634, 

.088) 

( 

.638, 

.769, 

-.032) 


The direction cosines for and are related to the direction 
cosines for and H 2 by a 180*^ rotation about the b axis. 

The smallest of the principal values for the 
coupling tensors reflected the largest uncertainties. These 
quantities represent the Interaction along the C-H bonds. 

Only ten of the thirteen parameters were determined 
with assurance that the parameters were actually Independent. Using 
all 3 principal values of the g tensor as parameters gave erratic 
results. In evaluating the Independence of only 2 g values as para¬ 
meters a dlagonallzatlon was performed on the D'D matrix found In 
the least squares calculation which Included the 2 g values as 
parameters. The large spread of the eigenvalues which resulted 
was Indicative of questionable Independence. This could have been 
avoided by choosing numbers of the same magnitude as the principal 
values of the coupling tensors and representing the g values In 
terms of these more appropriate numbers. However, a run was made 
using 12 parameters and the principal values of the g tensor became 
g - 2.00242 + .00004 

g - 2.00353 + .00006 

*yy 

g " 2.00420 was fixed and not varied. 
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TABLE II. Parameters obtained from least squares fit. 


PARAMETERS 


1 

2 

3 

4 

5 

Twist 


- 4.79 

- 5.26 

- 5.25 

- 6.30 

- 5.18 

Wag 


- 4.79 

- 5.00 

- 5.00 

- 5.35 

- 5.14 

Flap 


- 9.86 

- 9.88- 

— 9.88 

- 9.64 

- 9.95 

HCH 


116.77 

116.90 

116.90 

116.72 

116.79 

PRINCIPAL 

VALUES 







Tensor 


-26.93 

-26;80 

-26.80 

' -26.71 

-26.17 



-59.48 

-59.60 

' -59.59 

-59.80 

-59.29 



-92.24 

-92.14 

-92.15 

-91.74 

-92.52 

Tensor 


-25.75 

-26.42 

-26.41 

-28.55 

-26.17 



-60.06 

-58.96 

-58.98 

-59.73 

-59.29 



-92.37 

-92.78 

-92.78 

-92.68 

-92.52 

g Tensor 



2.00243 







2.00354 




STANDARD 

DEVIATION 









.417 

.335 

.335 

.426 

.339 


Parameters obtained with 

1. principal g values fixed at 2.0020, 2.0034, 2.0042 using 166 
points for the least squares fit. 

2. only 1 principal g value fixed; g « 2.0042; used 166 points. 

z z 

3. principal g values fixed at 2.00242, 2.00353, 2.00420; used 
166 points. 

4. Same as 3 except 244 points used for least squares fit. 

5. assumption that coupling tensors are same; principal g values 
fixed at 2.00242, 2.00353, 2.00420; used 166 points. 




























These values were then used in a program varying 


only 10 parameters to obtain the best fit possible. 

Due to different sets of data used and the number 
of parameters used for the least squares fit, the variation in 
the parameters was explored for these situations. The results 
are tabulated in Table II. One set of parameters was calculated 
using the assumption that both oc.-proton coupling tensors were 
the same. Table II also shows the use of all 6 principal values 
of the coupling tensor. The comparison of these parameters shows 
the validity of assuming that the coupling tensors are the same. 

The data taken for the 166 and 244 point programs was from the 
same spectra but the lowest number of points represents a refine¬ 
ment in this data. Differences were taken only between distinct 
and assignable lines. Overall the parameters are uniquely deter¬ 
mined within experimental error. The standard deviation for each 
set of parameters obtained Is Included in Table II. 

4.2 Spectrum of CH 2 CO 2 radical 

4.2.1 Observation of the "g-shift" 

The basic spectrum for interaction of an electron 
with two inequivalent protons is a quartet with equal intensities (28). 
To illustrate the "g-shift", the spectra for 2 field orientations 
are shown in Figure 10. Under these spectra is plotted the stick 
spectrum obtained from the solution of the spin Hamiltonian in 
PREDICT. The upper spectrum was taken with the magnetic field 
perpendicular to the y axis (or b axis). Since the b axis was 
established as the two fold axis in the crystal the two magnetic 
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10. ESR spectra of the ‘CH^CO* radical at 133°K. Magnetic field parallel to b axis (top) and 
perpendicular to z axis in the x-y plane (bottom). Predicted stick spectrum plotted for 
comparison. 
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sites should be equivalent and the observed quartet agrees with 
this prediction (12). The coupling constants are the same for 
both sites. The lower spectrum was taken with the magnetic 
field perpendicular to the z axis and the effects of shifting 
each site due to existing anisotropies are apparent. All of 
the spectra taken perpendicular to the y axis exhibited the same 
characteristic of equivalent sites. 

4.2.2 Symmetry of Inequivalent Sites 

The stereographic net can again be used to 
demonstrate the symmetry present in a single crystal of 

Referring to the Wulff net in Figure 6 it 
has already been mentioned that there are 79 different field 
orientations possible without duplicating spectra. If the 
right quadrant of the net is folded onto the left quadrant the 
spectra for these field directions will be the same but the 
sites will be interchanged. This is a result of the two fold 
symmetry about the b axis. And by the symmetry of a itereographie 
net the bottom half of the sphere is a replica of the top half, 

For a^J^j orientation the direction cosines were found to be 
(-.707,-.5, .5) by the following transformations 

axtp P«C. 

c sin 0 cos 0 

b -sin 0 sin 0 

a*i -cos 0 
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Figuro 11, Comparison of ESR spectra of the ‘GHjCO' radical at room (Cop) and Liquid nitrogen (bottom) temperature 
with magnotic field parallel to x axis, 

1 
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These equations vere direct results of the definitions of 9 and 

V 

0 in the mounting of the crystal In section 3.4. The spectra 
taken with field direction cosines of (-.707.5,.5) were the same 
as the spectra taken with the following direction cosines for the 
field: 

(+.707,+.5,-.5) 

(-.707,+.5,+.5) 

(+.707,-.5,-.5) 

When only the sign of the direction cosine for the b axis was 
changed the sites were Interchanged, but the spectra remained 
the same. This concept was particularly useful In manipulating 
the data. 

4.2.3 Comparison of Room Temperature Spectra with 
Low Temperature Spectra 

At room temperature the coupling constants for 
each proton were found to be the same with the magnetic field 
along the a+ axis. On cooling to low temperatures the radical 
encountered a barrier to Internal rotation about the C-C bond, 
and the center line split Into 2 lines (Figure 11). The room 
temperature spectra for the field along the b axis also showed 
a triplet which split into a quartet on lowering the temperature 
to 133^K. The sum of the coupling constants for the 2 protons 
should be the same at low and high temperatures but these sums 
taken from the two spectra were not the same. 
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TABLE III. 


Changes in the sum of the coupling constants 
(a^^ + a^^ ) for the methylene protons between 

12 ^ 
room temperature and 133 K, 


Field a + 

Orientation 1 2 

298°K 


"l “2 
133°k 


Change 

in 



O 00 

45.3 

47.5 

2.2 

ft’) 

42.0 

44.1 

2.1 

(?) 

38.9 

41.7 

2.8 

{“) 

36.8 

40.6 

3.8 

(?) 

36.6 

40.7 

4.1 

o o 

44.5 

47.7 

3.2 

/'90\ 

40.0 

42.6 

2.6 


48.5 

51.1 

2.6 

/90\ 

34.8 

40.6 

5.8 

[usj 

49.2 

52.7 

3.5 

(90\ 

36.9 

38.7 

1.8 

[is) 

44.2 

46.0 

1.8 

A5\ 

44.7 

45.4 

.7 

l3o; 

51.8 

55.9 

4.1 

/60\ 

41.3' 

42.6 

1.3 

[so) 

51.8 

55.5 

3.7 

(75) 

38.1 

40.8 

2.7 

(so) 

49.9 

53.4 

3.5 

/45\ 

41.1 

42.6 

1.5 

[eoj 

53.7 

56.5 

2.8 

*For each field 

orientation 

(a^ + a ) is given 

“Hi H 2 

for both sites 

the crystal. 
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Field $ 

Orientation 

K 

298°K 

*“i ““2 

1:^3 °K 

Change 

in 

feo) 

38.8 


41.7 

2.9 

(,60/ 

53.7 


56.8 

3.1 


37.3 


41.0 

3.7 

(^60/ 

52.1 


55.4 

3.3 

/60\ 

39.7 


41.9 

2.2 

(.75/ 

51.3 


55.3 

4.0 


39.5 


41.5 

2.0 

(. 75 / 

50.3 


53.5 

3.2 
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Work has been done on the isotropic proton 
hyperfine interactions in the ESR of fT-electron free radicals 
and it was concluded that out-of-plane vibrational motions 
should make a small contribution to the splitting as a result 
of the direct overlap between the proton and the unpaired 
electron in the TT-orbital (28) (1). A contribution such as 
this would decrease with increasing temperature because of the 
relationship between the proton and electron spin with respect 
to out-of-plane vibrations. As out-of-plane vibrations Increase, 
the sum of the coupling constants would then decrease. A 
theoretical analysis of the effect of vibrations on the methyl 
radical has been carried out by Schrader and Karplus (19). 
Zlochower, Miller, and Fraenkel applied this theory to the 
methyl radical and found a small dependence of hyperfine 
splitting with temperature (34) on the order of .04 gauss for 
a 20°C change in temperature in a 50°C range extending above 
and below room temperature. In this Investigation a series of 
spectra taken at room temperature were compared with those at 
133°K, and a change in hyperfine splitting of approximately 3 
gauss was found. Some of these differences are tabulated in 
Table III for various field orientations. 

4.2.4 Comparison of Calculated and Observed 2nd 

Order Lines 

Suspected 2nd order lines were readily observed 
in spectra taken with the magnetic field perpendicular to the b 
and c axes. These lines were more intense with the field perpenr 
dlcular to the b axis. Other directions of the magnetic field 
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Fi^jjrr: 12. Experimental (E) and predicted (P) spectra of the *0112002 showing 2nd order lines. 
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resulted in such an overlap of the spectra from the two sites 
that the 2nd order lines were either buried by more intense lines 
or distorted by shoulder effects. Figure 12 shows experimental 
spectra in which 2nd order lines are observed. The spectra predicted 
by PREDICT were included below each experimental spectrum. The 
predicted 2nd order lines are not as intense and the experimental 
plot shows lines not predicted. Strangely enough, PREDICT computes 
practically no 2nd order lines for the octant in the 

stereographic net while for nearly every orientation in the (+,-,+) 
octant, 2nd order lines are predicted with intensities ranging from 
,1 to .32 (relative to 1.0). The 45° bisector of this octant 
contains the orientations with the most Intense 2nd order lines 
that were predicted, 

4.2.5 Satellites 

A D^O substitution for H^O was made in the zinc 
acetate in an effort to explain the appearance of satellite lines 
in the spectra (Figure 13). In an attempt to study similar satellites 
in irradiated sodium acetate trlhydrate, Rogers and Klspert substituted 
D 2 O for H 2 O in the single crystal and observed a 6 fold reduction in 
satellite intensities (17). 

The ESR spectra of irradiated Zn(CH2C02)2‘D20 
crystals were taken at 133°K for selected orientations of the magnetic 
field (Figure 14). The crystals with D 2 O were analyzed by NMR. A 
sample with only H 2 O present resulted in an integrated peak of 88 mm 
and the sample with D 2 O showed an integrated peak of 38 mm which 
included some OH fragments. 
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Figure 13. ESR spectra of the radical. |90j field direction at room temperature. Other spectra at 133°K. 
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ESR spectra of .Cll^CO* radical from irradiated single crystals of Zn(Cll^CO^> 2 '^>2^ taken at room 
temperature (top) and 133 “k (bottom 2 spectra). 
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The spectra of the D^O species were found to be 
identical to those without the D^O substitution. One set of 
satellite lines are spaced approximately 5 gauss on either side 
of the main lines in both species, A second set of satellite 
lines are observed approximately 21 gauss from the main lines 
at high gain. The more intense satellite lines were found to 
increase in intensity at higher power but did not saturate. They 
appeared more prominent in spectra taken with the magnetic field 
positioned in the plane perpendicular to the two-fold axis. At 
these positions of the magnetic field the two sites were beginning 
to merge giving only one pattern for the two sites. Thus the 
satellite lines were less subject to shoulder effects and to 
dominance by the more intense lines of the quartet. 

The satellites in sodium acetate trihydrate were 
finally attributed to "spin-flip" transitions (17), which arose 
when there was a simultaneous change in the spin state of the odd 
electron on the methyl radical and of neighboring protons. Analogous 
results were a possibility in this investigation with the odd electron 
located on the oc-carbon of the methylene group and the spacing of 
the satellite lines occurring at a value close to the proton resonance 
frequency (5.0 to 5.8 gauss) in the magnetic field employed. Since 
no change in the spectra was observed using the D 2 O substitution the 
neighboring protons that could be responsible for the spin flips 
would have to be located in a species oriented near the radical. 

The protons of undamaged CH^CO’ were considered likely possibilities, 
but the most plausible theory in this investigation would be based on 
a spin flip of one proton concurrent with a spin flip of the unpaired 
electron. 
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Evidence o£ a single spin flip has also been 
observed in the spectra of gamma*irradiated single crystals of 
hydroxylated organic compounds (16) by Rao, et.al., using the 
theory of Trammel, Zeldes , and Livingston (27). This theory 
vas applied in this investigation with a proton placed near 
the methylene group in the * 011200 " radical but the resulting 
predictions were not altogether conclusive. The problems with 
shouldering effects, overlapping lines, and 2nd order lines were 
partially responsible for the unsuccessful analysis. And of 
course another species could be present making the assignment 
even more difficult, 

4,3 Spectrum of Irradiated Zn(CH2DC02)2*2H20 

Spectra were taken at 133*^K for selected field 
orientations since this investigation was primarily concerned 
with the Isotope effect of a proton abstraction to form the 
'CHDCO^ radical, A spectrum for one field orientation is shown 
in Figure 15, By orienting the magnetic field perpendicular to 
the twO"fold axis the g anisotropy was eliminated. This gave 
a spectrum that was easily comparable to a hyperfine pattern. 

A computer program was used to simulate the hyperfine pattern 
resulting from the substituted deuterium (24). The spectrum 
predicted by this program is also shown in Figure 15 to indicate 
the agreement obtained. The coupling constants taken from main 
line separations in the experimental spectrum compared favorably 
with those from spectra of the .CH 2 CO 2 radical for the same field 
orientation. 


67 








-1 

10 G^U&S 


ESR spettriim oT CHDCi)“ rat! u:a 1 with iiiagntdtic L'iGld per pend Lt'ii lar to b axis. Hot tern spectrmii prc'dict 
assuming a 2:L ratio of CIlDCO^ : CIl^CO'. 
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The hyperfine structure was predicted using the two 
radicals CHDCO^ and CH^iCO- but realizing that CHDCO- radicals 
would exist as conforraers. The coupling constant for deuterium 
was known to be .153 times that of hydrogen and since a deuteron 
may split each proton line into a triplet, the overlap of the 
lines in the center of the spectrum was to be expected. Also, 
this overlap was hard to interpret because of the presence of 
the CO^ radical in the center portion of the experimental spectra. 
However, the quartet frcxn the -CH^CO^ radical was easily identified 
and this formed a good basis of comparison between the predicted 
and experimental spectra. 

In considering an isotope effect two mechanisms 
involving the formation of •CHDCOg and •CH 2 CO- radicals were 
postulated. When a methyl radical (created by radiation damage) 
is located in a position close to an undamaged acetate ion, an 
isotope effect can exist; a C-H bond is more easily broken than 
a C-D bond so that the ensuing reaction 

iHiD +CHjDC0' = CHjD+'CHDC0^ 

results in formation of the •CHDCO^ radical. For a non-radical 
type reaction involving hydrogen and deuterium the maximum 
isotope effect has been reported as 6.9:1 (32). But with the 
methyl radical taking part in the reaction a large isotope effect 
is expected (31). If this is the case, the breaking of a C-H 
bond can be termed the barrier to formation of the radical. But 
the methyl radical may also migrate through the crystal until it 
approaches an acetate ion and randomly selects a hydrogen or 
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deuterium. Thus a statistical basis Is created when there Is an 


equal probability of selection-. The barrier to formation of the 
radical is then the migration of the methyl radical through the 
crystal. If the resulting abstractions are truly statistical 
two conformers of ‘CHDCO^ will be formed for each *0112002 radical. 

To show this statistical effect in the predicted hyperfine patterns, 
the Intensities for each radical were adjusted on a quantitative 
basis. The quartet was assigned a relative intensity of 1 and 
the two conformers whose pattern consisted of 6 lines each were 
assigned a value of .67. Thus the sum of the intensities for 
each radical was 4 and the statistical ratio of 2:1 (OHDCO^ over 
•0H2002) was maintained. The agreement obtained between predicted 
and experimental hyperfine patterns was excellent but the statistical 
model used was based on 100% deuteratlon. An analysis using ESR 
Intensities of the 'OH^ radical vs the *01120 radical from spectra 
of irradiated Zn (0H2D002)2'2H20 at 113°K showed only 79+3% 
deuteratlon of the acetate ion. If there were only 79% OH 2 DOO 2 
and 21% OH^CO^ present initially, the resulting ratio of OHDOO^: 
OH 2 CO 2 would be 1.1:1 assuming the abstraction of a proton to be 
entirely satlstical. But the experimental spectra still agreed 
with the predicted hyperfine pattern which assumed the 2:1 ratio. 

For a general approach: 





a 


b 


c 




1 -a 
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where = isotope effect 

a * b + c 

But experimentally 

r = - Y^-fc this investigation) 

Solution: 

b- - >r _ 

2e: 

If r ■» 2,0 ±. . 3 

a = .I*) i. .^>3 +Kcn ^ - 3.0 ±. 10 

2c 

So there is a significant Isotope effect and the barrier to 
formation of the •CH2C0^ radical is the breaking of the C-H 
bond on the acetate ion. The methyl radical probably resides 
in a position near the acetate ion from which it abstracts a 
proton. 
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5. Molecular Orbital Calculation for the CH 2 CO 2 Radical 

In working with the *0112002 radical the spin density of 
the electron was assumed to be concentrated on the o(.-carbon. 

The results obtained in section 4 showed this to be correct. 

To verify this a molecular orbital calculation based on the 
CNDO approximation was made (25)(15). The spin density on the 
ot-carbon-for the planar radical was found to be .69 with each 
of the oxygens having a density of .15. The calculation for 
the perpendicular form of the radical yielded a spin density 
of .82 on the o<.>rcarbon. The value of .69 for the planar 
configuration necessitated a check on the spin density for 
the o(.-carbon of the "CH(C 00 H )2 radical because the principal 
values of the coupling tensor for this radical had been used 
as initial -parameters in this investigation. The spin density 
on the ot-carbon of the ’ 011 ( 00011)2 radical was calculated to 
be .66 and in good agreement with the planar configuration 
of the ’0112002 radical. 

The barrier to Internal rotation of the methylene protons 
in the •OH 2 OO 2 radical was observed by Valenti and other investi¬ 
gators previously. The molecular orbital calculations on the 
planar and perpendicular models showed the planar form to be 
more stable. This degree of stability represents the barrier 
to internal rotation which turned out to be 11.8 kcal. The 
experimental value for the barrier has been calculated (using 
two different methods) to be 5.5 kcal (26). The difference 
between these values is somewhat attributed to the environment 
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In the crystal, but the errors inherent In such a calculation 


must be considered also. Nevertheless, a large zinc 
coordinating four formula units of the acetate In the 
which makes the lattice a unique home for the *CH 2 C 0 ^ 
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6. Conclusions 


A positive Identification of the *CH2C0^ radical In an 
Irradiated single crystal of zinc acetate dlhydrate has been 
made. The orientation of the radical In the crystal has been 
determined-with an HCH angle-of 116.9° +"4 3°. The radical's 
configuration Is nearly planar,- a result In'agreement with 
other Investigations- <2)' (22). The''principal values of the 
A and g tensors were found to'be of the' same magnitude as 
tensors of -other similar radicals In radiation damage studies. 

As In other TT -electron radicals the spin density was found 
primarily on the ec-carbon (20). A molecular orbital 
calculation supported this result and from the calculation 
the planar -form of the radical was observed to be the stable 
configuration. 

It has been observed In the-comparison'of spectra of 
Irradiated-single crystals of Zn(CH2DC02)'2’2H20 with predicted 
hyperflne patterns that a significant Isotope effect exists. 
Proton abstractions from acetate*ions by neighboring methyl 
radicals are considered the primary means for the formation 
of radicals. The barrier to formation of the *0112002 

radical Is then the breaking of the'O-H bond on the undamaged 
acetate Ion. 

A temperature dependence-of the coupling constants for both 
methylene protons has been observed In spectra at 298*^K and 133*^K. 
The magnitude of the coupling constants Increases with decreasing 
temperature and at a rate approximately 10 times larger than 
the rate found by Zlochower and co-workers In their study of 
the methyl radical (34). 
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The satellites appear"to'be "spin-flip"'transitions which 
arise from a change In the'spin'state of a proton on an undamaged 
acetate Ion concurrent with a change In spin state pf the electron 
on the o<.-carbon, but a further study of satellite line behavior 
Is warranted. The D^O substitution had no effect on the satellite 
lines; thus the water molecules are far enough removed from the 
unpaired electron to limit Interactions and It Is concluded that 
no dlpole-dlpole Interaction takes place between the water and 
methylene protons of the *CH 2 C 0 ^ radical. 

The 2nd order lines predicted by OPTZM are still not clearly 
assigned In the experimental spectra due to the confusion with 
possible spin flip satellites. There Is some evidence from 
spectra taken perpendicular to the b axis that the Intensities 
of the satellites are Increased due to superposition with 2nd 
order lines. 
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APPENDIX I 


A computer program for calculating transition frequencies 
and intensities for radicals oriented in a crystal lattice. 

The program is written in Fortran IV for use with the IBM 360 
computer. 
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VO 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM PREDICT 


NN IS 
IFRPT 


A. NUCLEI TO BE INCLUDED IN THE CALCULATIDJ^ 

INITIAL SET OF DIRECTION COSINES, THEN NOT ZERO FOP EACH 

nP^fuI^SArMcWr^lTcr'r.^^ CHANGING THE DIRECTION COSIjJeS^ 

Ur iHt nAbNcTlC FIELD 

^ NUCLEAr/^ZEEMAN°FREOUENCy'^^‘'’'‘“^*^ 'MEASURED IN MEGACYCLES OF THE 

COSINES MATRIX. ELEMENT I,J,K IS EOR NUCLEUS I, 

A PRINCIPAL AXIS WITH THE KTH CRYSTAL AXIS ’ 

PRINCIPAL VALUE OF THE COUPLING TENSOR FOR ATOM I 
UH(n IS THE I-TH COMPONENT OF THE UNIT VECTOR OOINTING ALONG THF 
DIRECTION OE THE MAGNETIC FIELD 

ARRAY OF TRANSITION FREQUENCIES CALCULATEO 
EINT IS THE CORRESPONDING ARRAY OF INTENSITIES 
IC IS THE NUMBER OF TRANSITIONS CALCULATED 


5 

10 

15 

16 
20 

25 


30 


31 

40 


50 


5,20 
1000,10 
NN 

((DCOSd, J,K) 

NN 

(A(I,J>, J = 1,3) 
(UH(I), I = 1,3) 


K = 1,3), J = 1,3) 


DIMENSION DCOS(10,3,3), 

REAO(5,100) NN,IFRPT,H 
IF(IFRPT) 5,“ 

IF(NN) 1000, 

DO 15 I = 1, 

REA0{5,101) 

DO 16 I = 1, 

READ(5,102) 

REA0(5,103) 

IFIIFRPT) 25,25,40 
WRITE (6,111) H 
WRITE(6,104) 

DO 30 I = 1,NN 
WRITE (6,105) ((DC0S(I,J,K), 
WRITE(6,106) 

00 31 I = 1,NN 

WRITE(6,107) (A(I,J), J = 1,3) 

WRITE(6,108) (UH(I), I = 1,3) 

CALL SSESR(NN,DC0S,A,UH,H,0E,FINT,IC) 


A(10,3), UH(3), OE(IOOO), FINT(IOOO) 


K = 1,3), J = 1, 


1000 

100 


WRITE(6,109) (DE(I), I 

WRITE(6,110) (FINT(I), 

DO 50 I = 1,IC 

DE(I) = DE(n/2.80246 

WRITE(6,112) (DE(I), I 

GO TO 1 

STOP 

FORMAT! 12,13,Flo.5) 


= 1,IC) 

I = 1,IC) 


= 1,IC) 




oono 


101 

102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 


FORMAT (9F7.5) 

FORMAT OFT,3) 

FORMAT OFT.5) 

FORMAT (IX,//* DIRECTION COSINES'//) 

FORMAT (1X,9F12.7) 

FORMAT (lX,//» COUPLING CONSTANTS'//) 

FORMAT (1X,3F10.5) 

FORMAT dXi/////' DIRECTION COSINES FOR H',//(3F12,7)//) 
FORMAT (IX,' TRANSITION FREQUENCIES',/(lOFlI.5)I 
FORMAT (IX,' TRANSITION INTENSITIES'./(lOFl1.5)) 

FORMAT ((IHl,'NUCLEAR ZEEMAN FREQUENCY = 'F11.5)//) 

FORMAT (//' TRANSITIONS, IN GAUSS FROM CENTER, ASSUMING G 
l',/(10F11.5) ) 

END 


2 . 002 ^ 


00 

o 


subroutine SSESR(NN,OCOS,A,UH,H,DE,FINT,IC) 

C NN IS THE NUMBER OF NUCLEI WITH SPIN = .5 

C DCOS IS THE DIRECTION COSINE MATRIX, ELEMENT I,J,K RELATES THE COMPONENT 
C FOR ATOM I OF THE J-TH PRINCIPAL AXIS ALONG THE K-TH CRYSTAL AXIS 
C A(I,J) IS THE J-TH PRINCIPAL VALUE FOR THE I-TH ATOM 
C UH(I) IS THE I-TH COMPONENT OF THE UNIT VECTOR POINTING ALONG THE 
C DIRECTION OF THE MAGNETIC FIELD 

C OE IS THE ARRAY OF TRANSITION FREQUENCIES CALCULATED 
C PINT IS THE CORRESPONDING ARRAY OF INTENSITIES 
C IC IS THE NUMBER OF TRANSITIONS CALCULATED 

DIMENSION DCOS( 10,3,3), A(10,3), UH(3), DEdOOO), FINTdOOO), 

1HV(10,3), HVP(10,2,3), CEdO), HC0S(10,3), ABH(10,2), S(2), 10(20) 

S(1) = 1.0 
S(2) = -1.0 
DO 10 I = 1,NN 
DO 10 J = 1,3 
HCOSd,J) = 0.0 
DO 10 K = 1,3 

10 HCOS(I,J) = HCOS(I,J)*UH(K»*OCOS(I,J,K» 

WRITE(6,100) ((HCOS(I,J), J = 1,3), I = 1,NN) 

DO 15 I = 1,NN 
DO 15 J = 1,3 
15 HVCI,JI = 0.0 
DO 20 I = 1,NN 
DO 20 J * 1,3 
CAS = HCOSCI ,J)^Ad,J) 

DO 20 K = 1,3 





2C 


30 


35 

40 


00 


45 

50 


55 


60 


64 

65 


66 

70 

100 


HV(I,K) = HVd.K) ♦ CAS*0C0Sf 
WRITEt6,100l ({HVfI,J}, J = 1*31, I = l,NNI 
00 30 I = IfNM 
00 30 J ^ 1,2 
DO 30 K = 1,3 

HVP(I,J,K) = H4UH(K)-0.54S(J)*HV(I,K) 

«RITEC6,100> f((HVP(, K = 1,31, J = 1,21, I 
00 40 I » 1,NN 
00 40 J * 1,2 
ABH(I,J) = 0.0 
00 35 K = 1,3 

ABH(I,J} = ABH(I,J) *■ HVPCI,J,K)**2 
ABH(I,J) = SORT(ABH(I,jn 

NRITE(6,100) ((ABH(I,J), J = 1,21, I » 1,NN) 


= 1,NN 
0.0 
= lt3 
CE(I) * 


HVPlI,l,JI*HVPtI,2,JI 


v,c«tf ••• nvri I , i.,u I »n*ri 1 , . 

CE(Il/fABH(I,ll4ABH(I,21} 
100) (CEd), I = 1,NN) 


DO 50 I 
CEIl) = 

00 45 J 
CE(I) = 

CE(1) = 

HRITEC6,iuv 
NN2 - 24NN 
NNTT = 2*4C2*NN) 

HNNT = 244NN 
00 55 I = 1,NN2 
10(1) = 1 
IC = 0 

00 70 I = 1,NNTT 
IC = IC41 
OE(IC) = 0.0 
FINT(IC) = 1.0 
00 60 J = 1,NN 
J1 == IQ(J) 

J2 = IO(NN+J) 

FINT(IC) = FINKIC)*(1.0-S(J1)4SCJ2)4CE(J)) 

OE(IC) = OE(IC) ♦ S( J1)4ABH( J,l) ♦ S( J2)4cABH( J,2) 
OE(IC) = 0.5*0E(IC) 

FINT(IC) = FINTCIO/HNNT 
J = 1 

IF(I0(J)>2) 66,65,66 

I0(J) = 1 

J * J+1 

GO TO 64 

I0(J) = 2 

CONTINUE 

FORMAT (IH ,10F11.5) 








uou 


RETURN 

END 


AN 

EXAMPLE 

OF THE 

DATA 

IS AS 

FOLLOWS 

2 

36.46 





-631 

299 

-715 

191 

955 

227 

968 

-141 

-213 

191 

955 

227 

290 

610 

910 




290 

610 

910 






1.0 




2 1 

36.46 





70711 
2 1 

36.46 

70711 




966 

2 1 

36.46 

259 




86602 

-50000 




2 1 

36.46 





70711 
2 1 

36.46 

-70711 





1.0 






-753 

170 


-Oil 

-265 


1 13.8 

21183 76829 

30 60 

1 

1 1 13.8 

1 

1 1 13.8 

1 1 13.8 

70711 70711 

1 1 13.8 

70711 
1 1 13.8 

70711 
1 1 13.8 

34202 93967 

1 1 13.8 

93967 34202 

1 1 13.8 

34202 
1 1 13.8 

93967 
1 1 13.8 

34202 


60402 -97091 
90 


09484 21985 11162 -63303 


70711 

70711 


93967 

34202 

93967 


659 

949 


76604 





APPENDIX II 


A computer program which-performs a least squares fit 
on selected data from predicted and observed spectra for a 
radical oriented in crystal'space; Differences between predicted 
frequencies are compared with the differences between observed 
frequencies and the least squares fit-'results in the determin¬ 
ation of certain selected parameters; The program is written in 
Fortran IV for use with the' IBM 360 computer. 
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C PROGRAM OPTIM 

c 

C THIS PROGRAM WILL CALCULATE ESR TRANSITION FREOS AND INTENS FOR AN 
C ORIENTED RADICAL IN A DEFINED CRYSTAL SYSTEM GIVEN THE PRINCIPAL 

C VALUES OF THE A AND G TENSORS. IT WILL PERFORM A LEAST SQUARES FIT 

C ON OIFFS TAKEN BETWEEN SELECTED FREQS FOR OBSERVED AND CALCULATED 
C SPECTRA. THE PARAMETERS ARE CORRECTLY DETERMINED BY THE LEAST 
C SQUARES FIT. 

C 

IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION ACOS(2» 3,3)♦A(2,3),UH(3),DE(16),FINT(16),PAR(13), 

1 CC(3),HH(3),UHM(70,3),DE0B(250),IFS(70,13),R(250), 
2DT(13,250),DEL(13),PARN(13),DEl(250),DE2(250),B(2,3),DEN(13 ) , 
3IR(70),DECOS(2,3,3),GCOS(2,3,3),GECOS(2,3,3),G( 2,3),OEE(16), 
4DE4(250),DE3(250),0E5(250),DE6(250),G1(2,3),NFR(7C),FINTE(16) 
5,IFE(70,13),NFRE(70) 

6,ID(60),LGR(7C),IFM(70,13),IFN(70,13),IS(70) 

7,IT(70) 

REAL D(32),FIN(32),YI,ZI(60) 

REAL+8 ITITLE(12) 

READ IN - NN = NO. OF PROTONSINN MAX = 2) 

H - NUCLEAR ZEEMAN FREQ 

DEN = THE INCREMENTAL CHANGE IN THE PARAMETERS 
USED TO OBTAIN THE NUMERICAL DERIVATIVES. 

R£A0(5,100) NN,H,OEN 

READ IN - CC £ HH = DIRECTION COSINES FOR VECTORS ON WHICH 

INITIAL ORIENTATION OF RADICAL IS CONSTRUCTED. 

READ(5,114) CC,HH 

READ IN - PAR = INITIAL PARAMETERS 

READ(5,113)(PARI I),I = 1,13) 

WRITE(6,111) H 

THE NEXT TWO ARRAYS ARE FOR SUBROUTINE SPEC 

READ IN - ID - CONTROLS NO. OF SPECTRA DRAWN ON EACH GRAPH. 

ZI - CONTROLS THE SPACING OF THE SPECTRA ON EACH GRAPH. 

READ!5,126)(I0(I),1=1,60) 

READ(5,128) (ZKI) ,1 = 1,60) 




c 

c 

c 

c 

c 

c 

c 

c 

c 

NOW BEGIN ITERATION PROCEDURE 

READ IN - NFLP = NO* OF MAGNETIC FIELD DIRECTIONS 

NP = NO. OF PARAMETERS OR VARIABLES FOR LEAST SQUARES 

IP * INDEX WHICH IF = 0 READS IN NEW SET OF 

parameters at END OF DATA AND IF = 1 PROCEEDS 
WITH ORIGINAL SET READ IN ABOVE 

c 

c 

c 

c 

c 

c 

c 

c 

c 

1 READ(5,117» NFLP,NP,IP 

REAO(5>15M ITITLEf Ilf I»l»12) 

READ IN THE NEXT THREE ARRAYS TO TELL THE COMPUTER hOW MANY 

FREQS (BETWEEN WHICH DIFFS WILL BE TAKEN! HAVE BEEN SELECTED 

FROM EACH SPECTRUM 

NFR SITE 1 

NFRE SITE 2 

LGR SITE 1 £ SITE 2 

CD 

^ c 

c 
c 
c 
c 
c 
c 

REAOISfllTMNFRd 1,1=1, NFLP) 

READ(5,117)(NFRECI),I=1,NFLPI 

READ(5,117) (LGRd ),t=l,NFLP) 

READ IN THREE MORE ARRAYS TO TELL COMPUTER TO INCLUDE (1) 

OR EXCLUDE (0) A PARTICULAR SITE FOR FREQ SELECTIONS 

IT SITE 1 

IR SITE 2 

IS SITES 1 £ 2 

REA0(5,117)dT(I),I = l,NFLP) 

REA0(5,117)(IR(I),I=1,NFLP) 

REA0(5,117)dS(I),I=l,NFLP) 

NF2=C.O 

Jl=1.0 

DO 200 1=1.NFLP 

NFE=NFREII) 

LG=LGR(I) 

ooo 

NF=NFR(n 

NF1»IT( I )*(NF-ll*IRd )*(NFE-l )♦!£(! I’KLC 

NF2»NF2*NF1 

UHM = MAGNETIC FIELD VECTORS 

REA0(5,103) (UHMd,J), J=l,3) 
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READ IN ARRAYS IFSflFE.IFM, IFN TO INDICATE WHICH FREOS TO 
SELECT WITHIN EACH SITE. IFM AND IFN FREOS ARE PAIRED FOR 
OIFFS BETWEEN EACH SITE. 

IFS SITE 1 

IFE SITE 2 

IFM SITE I 

IFN SITE 2 

REA0(5,126) (IFSIItJ)tJ^ltNF) 

READ!5t1261(IFE(ItJ)tJ^ltNFE) 

REAO(5tl26)nFN(I»J)»J«l,LG) 

REAO(5f 126) (IFNn»J)»J=l,LG) 

READ IN experimental OIFFS TO COMPARE WITH CALC FREOS. 

READ!5*113) (OEOBU) tJ-JliNF2) 

WRITE(6*108)(UHM(I»J)•0=1*3) 

WRITE(6*113)(0E0B(J)*J=J1*NF2) 

WRITE(6*129) NF2 
WRITE(6.118) NF1*I 
200 J1=J1+NF1 

IF(IP) 2*2*3 

2 REA0(5*113)(PAR(I).I=1,13) 

WRITE(6*116)(PAR(N)*N=i*NP) 

3 CONTINUE 

WRITE(6*115)<PAR(I)*I=1*NP) 

START ITERATION FOR LEAST SQUARES FIT. 

DO 60 M=l,3 
WRITE(6*106) 

DO 16 J=l*3 

J4=J+A 

J7=J+7 

J10=J+10 

Ad* J) = PAR(J4) 

A(2*J)=PAR(J7) 

16 G(1,J)=PAR(J10) 

DO 31 I = l.NN 

31 WRITE(6*107) (A(I*J)* J = 1*3) 

WRITE(6*124)(G(1*J)*J=1*3) 

ENTER TWSF TO UPDATE DIRECTION COSINES FOR EXISTING PARAMETERS. 





c 

c 

c 

c 

CALL TWSF(CCtHHtAC0S,0EC0S»GC0St6EC0S»PAR) 

WRITE(6f1041 

DO 30 I - 1,NN 

30 gglTE ((ACOS(I»J,K)» K = lt3)» J = 1,3) 

33 WRITE(67i65)((OECOS(I,J,K),Ks1,3),J»1,3) 

MRITE(6,125)( (GCOSIl. J,K),K-1,3 ),J:s1,3| 

MRITE(6,125)(fGEC0SCi,J,K),Ksl,3),Jsl,3) 

L=i 

JJ=1 

KKsl 

00 20 J=1,NFLP 

DO 35 K=l,3 

35 UHCK)=UHM(J,K) 

ENTER SSESR TWICE FOR PREDICTED FREOS AND INTENS. FOR BOTH 
SITES!MONOCLINIC CRYSTAL) 

00 

•sj 

CALL SSESR!NN,ACOS,GCOS,A,G,UH,H,OEiGEFfFINT,IC) 

CALL SSESR(NN,OECOS,GECOSfA,G,UH,H,DEE,GEFE,MNTE,IC) 

D0~36 1=1,IC 

0(IX)=0EE(I) 

FIN(IX)=FINTECI) 

IX=IX+1 

0(lX)=0E(I) 

FIN!IXI=FINT(I) 

36 IX=IX+1 

NF=NFR(J) 

MODCUR=IOCJ) 

YI=ZICJ) 

c 

c 

c 

c 

IF DESIRED ENTER SPEC ON SELECTED ITERATION CM) FOR DRAWING 
PREDICTED SPECTRA 

IFCM-4) 38,37,38 

37 CALL SPEC(ITITLE,D,FIN,IC,MODCUR,YI) 

38 CONTINUE 

c 

c 

c 

BEGIN FREQ SELECTION FROM EACH SITE 

DO 40 K=1,NF 

MM=IFSCJ,K) 

OEK L)=OE(MM) 





oooo 


C 


40 L=L+1 
NFE»NFRE(J) 

00 41 K=1,NFE 
MM=IFEfJ,K> 

DE3(JJI=0EE(MM| 

41 JJ=JJ+1 
LG»LGR(J) 

00 42 K»1,LG 
MM=IFM<J,K) 
KL=IFN(J,K) 
DE2(KK|=DE(MM| 
DE5(KK»=DEE(KLI 

42 KK=KK-H 
LGF=KK-LG 
K1*KK-1 
LME=JJ-NFE 
LE=JJ-1 
LNF=L-NF 
LM=L-1 
I1=LNF^1 
I2*LNE«-1 


TAKE DIFFERENCES FROM SELECTED FREQS FOR SELECTEO SITES AND 
STORE UNTIL CALCS USING ALL FIELD VECTORS ARE COMPLETE 

IFfIT(J)» 47,47,48 
48 CONTINUE 

00 410 N*I1.LM 

0E4CLL)=DABS(0E1(LNF|-DE1(N)I 


410 LL=LL+1 
47 CONTINUE 

IF(IR(J)| 46,46,45 

45 CONTINUE 

00 70 II=I2,LE 

0E4< LLI »OABS (0E3(LNE)-0E3 (III) 
70 LL=LL+1 

46 CONTINUE 
IFdSCJI) 20,20,44 

44 CONTINUE 

DO 71 L1»LGF,K1 
DE4(LL)=0ABS(DE2(L1)-DE5(L1)I 


71 LL=LL+1 
20 CONTINUE 
NFD=LL-1 
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CALCULATE R MATRIX WHICH IS CALC - OBS OIFFS 

DO 25 1=1,NFO 
25 R( I )=0E4( n-DEOB( I » 


SET UP DERIVATIVE MATRIX BY 
TO OBTAIN CHANGE IN FREQS. 
DETERMINING R MATRIX. 


VARYING EACH PARAMETER SEPERATELY 
AND PROCEDURE IS SAME AS IN 


00 

VO 


51 


61 


53 


54 


56 


DO 50 J=1,NP 
DO 51 K=l,13 
PARN(K)=PARIK) 

PARN( J» = PAR(JH-OEN( J» 

DO 61 K=l,3 
K4=K+4 
K7=K+7 
K10=K+10 

GKl ,KJ=PARN(K10) 

B(1,K)=PARN(K4» 

B(2,K)=PARN(K7) 

CALL TWSF(CC,HH,ACOS,DECOS,GCOS,GECOS,PARN) 
KK=1 


LL=1 

L=l 

JJ=1 

00 52 N=1,NFLP 
DO 53 K=l,3 
UH(K)=UHM(N,K) 

CALL SSESRINN,AC0S,GC0S,8,Gl,UH,H,OE,GEF,FINT,IC» 
CALL SSESRINN,DECOS,GECOS,B,G1,UH,H,DEE,GEFE,FI NT 
Nf=nfr(NI 

DO 54 K=1,NF 
MM=IFSIN,K) 

OEK L)=OE(MM) 


L=L+1 

NFE=NFREIN) 

DO 56 K=1,NFE 
MM=IFE(N,K) 

0E3(JJ)=DEE(MM» 

JJ=JJ+1 

LG=LGR(N) 

DO 57 K=1,LG 
MM=IFM(N,KJ 
KL=IFN(N,K) 
DE2(KKI=DE{MM) 


IC) 






r>oo oooo ooo 


o 


57 


216 


710 

215 

210 


220 

58 

224 


225 

52 


0E5(KK)»0EE(KL) 

KK=KK*1 

LNE=JJ-NFE 

LE=JJ-1 

LGF=KK-LG 

K1=KK-1 

LNF=L-NF 

LM»L-1 

I1=LNF*1 

I2=LNE+1 

IF(IT(N)) 215»215f216 

CONTINUE 

DO 710 K=I1.LH 

0E6(LL)=0ABS{DEI(LNF)-DElfKn 

LL=LL*1 

CONTINUE 

IF(IR(Nn 58,56t210 

CONTINUE 

00 220 II-I2.LE 

0E6(LL»=DABSCDE3CLNEI-0E3(111) 

LL=LL+1 

CONTINUE 

IFflSfNM 52t52t224 

CONTINUE 

00 225 L1=LGF,K1 

DE6lLL»=0ABS(0E2(Ll»-0E5ILin 

LL»LL+1 

CONTINUE 


NFO = NO. OF OIFFS TAKEN 


NFD=LL-1 
DO 50 K=1,NF0 

50 DT(JtK) = (DE6fK)*0E4(KH/DEN(J) 

THE DERIVATIVE MATRIX COTI HAS BEEN FORMED. ENTER BSTFT 
FOR LEAST SQUARES FIT 

CALL BSTFT(0T,R,0EL,NP,NF0l 
DO 55 K=1,NP 
55 PARCKI=PARCK|-OEL(K» 

60 WRITEf6,I16MPAR(N}fN=l»NP} 

PAR CONTAINS NEW PARAMETERS FOR NEXT ITERATION 




oooo 


GO TO I 

15 FQRMAT(6A8/3A8/3A8I 

100 FORMATn5,F10*5,13F5.0l 

101 FORMAT <9F7.5) 

102 FORMAT (3F7.3) 

103 FORMAT (3F7.5I 

105 FORMAT* ^IX^9F12*7|^'*'*°^ COSINES FOR THE "A** PRINCIPAL VALUES'//) 

106 FORMAT (lxI//» COUPLING CONSTANTS'//) 

107 FORMAT (1X,3F10.5) 

108 FORMAT (IX,/////' DIRECTION COSINES FOR H',//(3F12.7)//) 

109 FORMAT (IX,' TRANSITION FREOUENCIES•,//(lOFll.5)//) 

110 FORMAT (IX,« TRANSITION INTENSIT IES»,//(lOFll,5)//) 

111 FORMAT ((IHI,'NUCLEAR ZEEMAN FREQUENCY = *F11,5)//) 

112 FORMAT (//' TRANSITIONS, IN GAUSS FROM CENTER, ASSUMING G = 2.0023 
!',//(10Fll.51//) 

113 F0RMAT(8F10.5» 

114 F0RMAT(3F10.6/3F10.6) 

115 FORMATdX,//' REFERENCE PARAMETERS',//(lOFll. 6)//) 

117 FORMATlisis) * PARAMETERS',//( lIFll.5)//) 

118 FORMAT! iXj//'^ ^OIFFS PER FIELD ORIENTATION = *15,' NO. 

119 FORMAT! 


121 

122 

123 

124 

125 


127 

128 
129 


*15//) 
X,//' G 


FIELD ORIENT 


1F10.5// 

120 FORMAT! 

FORMAT! 

FORMAT! 

FORMAT! 
FORMATdX,//' 
FORMATdX,//' 
110F11.5)//) 

126 F0RMAT!24I3) 
FORMATdX,/' r 
FORMAT!24F3.0) 
F0RMAT!4X,/' TOTAL 
END 


EFFECTIVE = 'F10.5,3X,' 6 EFFECTIVE,E0IV.SITE = ' 


X*//(10Fll**i)//)'■*' " FRE0S',//!10F11. 5)//) 

xj//' OT(I*J): CHANGE IN PARAMETER PER CHANGE IN POINT',/) 
X,//' CALCULATED DIFFERENCES',//dOFl 1. 5)//) 

THE PRINCIPAL VALUES OF THE G TENSOR ARE'//!3F10.5)) 
DIRECTION COSINES FOR THE "G" PRINCIPAL VALUES'//! 

FOR THE EQUIVALENT SITE'//) 

POINTS OBSERVED = ',15,//) 


C 

C 

C 


SUBROUTINE SPEC!ITITLE,D,FINT,IC,MOOCUR,YI) 

SPEC DRAWS PREDICTED SPECTRA USING A LORENTZIAN DIST'N 
PLACES 4 SPECTRA ON A GRAPH. MAX OF 15 GRAPHS. 
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c 


DIMENSION F(900),0(321,FINT(32) 

1vW(900) 

REAL*8 ITITLE(12) 

W(1)=0.000 
DO 16 1=1,899 
IX=I+1 

16 W(IX)=W(I)+.100 
DO 71 1=1,900 
71 F(I)=0,0 
ICT=IC^2 
DO 120 1=1,ICT 
120 D( n=D( I)+37.5 

WI = LINEWIDTH 

WI=2.1 

T2=2/(1.73502*WI) 

DO 72 1=1,900 
DO 72 J=1,ICT 

F(I)=F( I)-2*(W(I)-D(J) )*T2*FINT(J)/(1+(T2+{W( I)-D( J) ) )«*2)*+2 

Uu liO 1—IfSOO 

W(I)=.877*W(I) 

110 F(I)=-.025+F(I)+YI 

^CALL^DRAW(900,F,W,MODCUR,0,LABEL,ITITLE,2.0E-02,6.,0,1,0,2,9, 

WRITE(6,80) LAST 

80 FORMAT!//* GRAPH OPERATION*/FIO.5) 

RETURN 

END 


C 

C 

C 

C 

C 

C 

C 

C 


SUBROUTINE BSTFT(DT,R,DEL,IR, IS) 


- OT TIMES D, 

OT IS THE TRANSPOSE OF THE DERIVATIVE MATRIX. 

TM^DAD AUCTCDC**^^^ PERFORMED GIVES BEST CHANGES TO BE MADE 
IN PARAMcTEKS« 

UNCERTAINTIES FOR EACH PARAMETER ARE CALC. 


IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION DT(13,250),8(250),DEL!13) 


15,0, 





ItOPIC13,l?),U(13.13),nXll3>,OPS(13,131 
PRINT IOh, hr,IS) 

PRINT 10*^ 

00 2C I = 1,IR 

PRINT 101, (DT(I,J), J = I,IS) 

20 PRINT 105 
PRINT 105 
PRINT 110 

PRINT 102, (R( I ) , I = 1,IS) 

PRINT 105 

CALL FIT( IR, IS,DT,R,DEL,DPI,U,DPS) 

RSUM = 0.0 
DO 30 I = 1,IS 
30 RSUM = RSUM+Rn)**2 
PRINT 106, (RSUM) 

AR=IS-1 

RSUM=RSUM/AR 

SD=DSQRT(RSUM) 

WRITE{6ill2) SO 
PRINT 105 
DO 35 I = 1,IS 
DO 35 J = 1,IR 

35 R(I) = R(I)-DT(J,I)*DEL(J) 

PRINT 111 

PRINT 102, (R{I), I = 1, IS) 

RSUM = 0.0 
00 40 I = 1,IS 
40 RSUM = RSUM4-R( I )**2 
PRINT 107, (RSUM) 

RSUM=RSUM/AR 

SD=DSORT(RSUM) 

WRITE(6,113) SO 
DO 45 1=1,IR 
DSX=RSUM»DPI(1,1) 

45 DX(I)=DSORT(DSX) 

PRINT 114 

PRINT 103, (DX(I),1=1,IR) 

DO 50 1=1,IR 
DO 50 J=1,IR 
50 U(J,I)=U(I,J) 

PRINT 105 
PRINT 108 

PRINT 103, (OEL(I), I = 1,IR) 

100 FORMAT (315) 

101 FORMAT (10F8.2) 









oooo 


102 FORMAT (10F8,4) 

103 FORMAT (IH , 5E20,8» 

104 FORMAT <//• NUMBER OF PARAMETERS * '12** NUMBER OF POINTS * *13, 

l//» 

105 FORMAT (IH » 

106 FORMAT (37H ORIGINAL SUM OF SQUARES OF ERRORS = tE15.6) 

107 FORMAT (32H NEH SUM OF SQUARES OF ERRORS = ,E15.8) 

108 FORMAT <41H ADD TO ORIGINAL PARAMETERS THE FOLLOWINGI 

109 FORMAT I* THE TRANSPOSE OF THE DERIVATIVE MATRIX IS*,//I 

110 FORMAT (//,* THE ORIGINAL ERRORS BETWEEN CALCULATED AND OBSERVED 
1ARE*,//I 

111 FORMAT (//,* THE ERRORS AFTER THE ITERATION ARE*,//! 

112 FORMAT!/,• ORIGINAL STD DEVIATION = •tElS.S} 

113 FORMAT!/,' NEW STD DEVIATION = •.EIB.BI 

114 FORMAT!//* THE UNCERTAINTIES IN THE PARAMETERS ARE*,//) 

RETURN 

END 


VO 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE FIT!IR,IS,DT,R,OEL,OPI,U,DPS) 

IR IS NUMBER OF VARIABLES 
IS IS NUMBER OF OBSERVED POINTS 

DT IS THE DERIVATIVE MATRIX, IR ROWS. IS COLUMNS 

R IS THE DIFF MATRIX, COLUMN, WITH IS ELEMENTS 

DEL IS THE VALUE TO BE ADDED TO EACH VARIABLE, A COLUMN MATRIX 

WITH IR ELEMENTS 

THIS PROGRAM SUPPLEMENTS BSTFT AND PERFORMS THE MATRIX 
MULTIPLICATIONS INDICATED IN BSTFT. 

IMPLICIT REAL*8!A-H,0-ZI 

DIMENSION 0T!13,250),DEL!13),R!250),DTR!13),0P!13,13),DPI!13,13) 
2.0PS!13,13),U!13,13) 

00 30 I = l,IR 
DTRil) - 0.0 
DO 30 J = 1,IS 

30 OTR!I) = DTR!I)+OT!I,J)*R!J) 

DO 40 I * l,IR 
DO 40 J = 1,IR 
0P!I,J) = 0.0 
DO 40 K = 1,IS 


40 0P!I,J) = DP! I,J)^OT!I,K))«'OT!J,K) 
DO 60 1=1,IR 





nooo oooo 


VO 


00 60 J=l,IR 
60 DPS<I,JI^DP(I.J» 

CALL HDIAG(DPS,U,0fIR,NR) 

, , ^ WRITE(6,112M0PS( I, n ,I = 1,IR» 

112 FORMATcIx,//' EIGEftvAtuEsJt//(l0E12.5)//) 

113 FORMATCIX,//* EIGENVECTORS».//(lOFl!•51//1 
CALL^GAysi3nR,1.00E-30,DP,6pi;iTEU» ' 

00 55 1=1,IR 

55 PRINT lll,(DPI(I,J),J=1,IR| 

111 FORMAtIiOEIuI) DERIVATIVE MATRIX IS'//» 

IFflTEST-l) 42,45,42 
42 PRINT 100 
45 00 50 I = 1,IR 
DELI II = 0.0 
DO 50 J = l,IR 

50 OELII) = DELII)+OPI(I,J|*OTR(J| 

100 FORMAT I16H MATRIX SINGULARI 

END 


SUBROUTINE HOIAGIH,U,IEGEN,N,NR I 

THIS ROUTINE DIAGONALIZES H AND DETERMINES THE ROW FrKFNVPrrnoc ii 
THIS VERSION OF HOIAG WORKS ON THE FULL MATRIX, NOT JUST THE UPPER HALF. 

IMPLICIT REAL*8IA-H,0-Z) 

DIMENSION HI13,13I, UI13,13I, XI13I, I0I13I 

M = 1 

FN = N 

DO 10 I = M,N 
DO 10 J = M,N 
10 UII.J) = 0.0 
00 12 I = M,N 
12 UII,II = 1.0 

15 NR = 0 
SUM = 0.0 

00 16 I = M,N 
DO 16 J = I,N 

16 SUM = SUM-t-HI I, J }>tc42 






o^ 


19 

20 
22 


23 

25 

39 

40 

41 


42 

43 


150 


SUM=0SQRT(SUM/(FN*FN>I 
SUM = SUM*l.00E-08 
N1 » N~l 
N2 = N-2 

IFCN2) 1000,40»22 
00 39 I * M,N2 
K = I+l 
XMAX 0*0 
DO 25 J = K,N 

lF(OABSfH(ItJ)l-XMAX) 25»25»23 
XMAX-OABS(H(I»J)) 

IQS = J 
CONTINUE 
IQ(I) s IQS 
Xdl - XMAX 
XINII s H(NltN) 

IQINII » N 

ZHAX » 0.0 

DO 43 I = M,N1 

IF(X(n-ZMAX) 43,43,42 

ZMAX s x(n 

IPIV * I 

CONTINUE 

JPIV = lOUPIV) 

IFIZMAX-SUMI 1000,1000,150 


NR s NR^l 
Z1 = HdPIVtIPIVI 
Z2 = HlJPIVtJPIVI 
Z3 = HdPIV,JPIV» 

Z4 * 4,0*Z3**2 
Z5 » (Z1-Z2I**2 
Z6 = (Z2-Z1»*Z3 
IF<0ABS(Z5/Z4»-10.0I 50,52,52 
50 SINSQ=0.5-0.5*0S0RT(Z5/(Z4+Z5II 
GO TO 55 
52 TC2 = Z4/Z5 

SINSO = 0*254TC24Cl,0-0*75*TC2»d.0-0,833333333333*TC2*(l,0-0.875* 
lTC2*d#0-0«9»TC2*(1.0-0#91666666667*TC2*( 1«0-0.92857142857*TC24 
2d.0-0.9375*TC2n dill 
55 COSSO = 1#0-SINS0 
COSINE*OSORT(COSSQ» 

SINE-OSQRT(SINSO) 

IF(Z6) 90,90,85 
85 SINE = -SINE 
90 SCH = 2.0*SINE*C0SINE*Z3 





|k 


SO 


100 


200 


201 


205 

210 

219 

220 


229 

230 

236 

238 


240 

245 


HIPIP = C0SSQ4Z1+SINS0*Z2*SCH 

HIPI = HdPIVyl) 

HJPI = HfJPIVfU 

HdPIVfl) = C0S1NE4HIPI^SINE*HJPI 
: 5f}gf;H{p.»cosiNE.HjP, 

H( I, JPIV) = H( JPIV.n 
H(IPIV,IPIVI = HIPIP 
h(JPIV,JPIVl = HJPJP 
HUPIV.JPIVI = 0.0 
HIJPIV,IPIVI = 0.0 
lOF = JPIV-IPIV 

ipm^I°iA¥Ii" ^p*v,jpiv,idf 

IF(IPM-N) 220»201,201 
XMAX =0.0 
00 210 I = IPM.N 

IF(OABS(H(INT,fn-XMAX» 210,210*205 
XMAX=OABS(Hf INTflH 

CONTINUE 
XIINT) = XMAX 
IQIINTI = IQS 
CONTINUE 

00^260jINT^= IPIV,JPIV»IDF 

IFtIPMI 260,260,229 
00 250 I = MtIPN 
IF(IO(II-INT) 240,230,240 

XMAX=DABS<H(I,INTII 
IQS = INT 
00 238 J = K,N 

IQS = J 
CONTINUE 
X(I) = XMAX 
IQCII = IQS 
GO TO 250 

IQCII = INT 
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250 

260 

500 


550 

555 

1000 


643 

645 


VO 

00 


650 

655 


CONTINUE 

CONTINUE 

IFCIEGEN) 41,500»41 
00 550 I = M.N 
UIPI » UUPIVtl) 

UJPI s U(JPIV»n 

U(IPIV,n = C0SINE*UIPH-SINE4UJPI 

UCJPIVtl) * -SINE*UIPI+COSINE*UJPI 

GO TO 41 

Nl = N-1 

00 655 J - 1»N1 

XNAX » HiJtJ) 

JQ = J 
Jl = J+l 
DO 645 I - J1»N 
IF(H(Itn-XMAXl 645t645,643 
XHAX s H(Itl) 

JQ » I 
CONTINUE 
TEMP = HIJ.JI 
HfJfJ) « H(JQ,JO> 

H(JQ.JO) » TEMP 
00 650 I M.N 
TEMP = U(J,n 
U(J«II > U(JQ.n 
UUQ.n * TEMP 
CONTINUE 
RETURN 
END 


SUBROUTINE GAUSS3(NtEPtAtX,KER) 
IMPLICIT REAL48(A-H,0-Z» 

THIS PROGRAM INVERTS A MATRIX. 

DIMENSION AU3,13)»X(13tl3) 

00 1 1=1,N 

00 1 J=1,N 

1 X(I,J)=0.0 
DO 2 K=l,N 

2 X(K,K»=l.O 
10 DO 34 L=1,N 




• « 


S w 




VO 

VO 


KP=0 

Z=0.0 

00 12 K=L,N 

,, IF(Z-OABSIAfK,L}n 11,12,12 

11 Z=DABS(A(K,Ln 

K P-sK 

12 CONTINUE 

,, IF{L-KPI13,20,20 

13 00 14 J-L,N 
Z:^A(L,J) 

A(L,J)»A(KP,J) 

14 A(KP,J|sZ 
00 15 J=1,N 
Z=X(L,JI 

XfL,J)=X(KP,J) 

15 XCKP,J)=Z 

20 IF(0ABS(AfL,Ln-EP} 50,50,30 

30 IF(L-NI 31,34,34 

31 LP1=L+1 

DO 36 K^LPl.N 
IF(A(K,L))32,36,32 

32 RATIO=ACK.LI/A(L,LI 
OQ 33 J=LP1,N 

33 A(K.JI-A(K,J)-RATI04A(L,J) 

DO 35 J=1,N 

35 X(K, JI=X(K,JI-RATI04'X(L,J) 

36 CONTINUE 

34 CONTINUE 

40 00 43 1=1,N 
II=N+1-I 

DO 43 J=1,N 
S=0,0 

IFni-NI41,43,43 

41 IIP1=II*1 

00 42 K=IIP1,N 

42 S=S*AUI,K»*X(K,J» 

43 X(II .J|vr(X(II,j)-S)/A{II,in 
KER=1 

RETURN 
50 KER=2 
RETURN 
END 


C 

C 

C 










100 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE SSESRINN»ACOS»GCOSfA,GtUH«H«OEtGEFtFINT»IC) 

THIS PROGRAM CALCS FREQS AND INTS* FOR ORIENTED RADICALS. 

ALSO CALC. EFFECTIVE G VALUE (GEFI AND COMPUTES THE CHEMICAL 
SHIFT OF EACH SITE RELATIVE TO EACH OTHER. 

NN IS THE NUMBER OF NUCLEI WITH SPIN = .5 

ACOS G GCOS ARE THE DIRECTION COSINE MATRICES FOR A AND G TENSORS. 
ELEMENT I.JfK RELATES THE COMPONENT FOR ATOM I OF THE J-TH PRINCIPAL 

AXIS along the k-th crystal axis 

AdtJI IS THE J-TH PRINCIPAL VALUE FOR THE I-TH ATOM 
UHin IS THE I-TH COMPONENT OF THE UNIT VECTOR POINTING ALONG THE 
DIRECTION OF THE MAGNETIC FIELD 
DE IS THE ARRAY OF TRANSITION FREQUENCIES CALCULATED 
PINT IS THE CORRESPONDING ARRAY OF INTENSITIES 
IC IS THE NUMBER OF TRANSITIONS CALCULATED 

IMPLICIT REAL*8(A-H.0-Z) 

DIMENSION ACOS(2»3,3)tA(2,3)»UH(3),0E(16),FINT(16). 

1HV(2,3)•HVP(2,2»3)tCE(10),HC0S(2t3),ABH(2»2)»S(2)»IQ(20)t 
2GCOS(2»3.3),GI2«3).PCOS(2,3) 

Sill ~ 1.0 
S(2)*-1.0 
DO 10 I»1.NN 
DO 10 J=1.3 
HCOS(ItJ) - 0.0 
PCOSdt J)»0.0 
DO 10 K^l»3 
PCOSd,J) = PCOS 
HCOS(I.JI^HCOS 
DO 15 i = 1*NN 
DO 15 J = 1,3 
HVd.J) s 0.0 
DO 20 I ^ 1,NN 
DO 20 J = 1,3 
CAS > HCOSf i.J)*Ad,J) 

DO 20 K = 1,3 

HVd,K) = HVd,K) ♦ CAS*ACOSd,J,K) 

I»1 

GEF=0.0 
DO 32 K=l,3 

GEF=6EF+Ck0S(I,K)*GCI,K) 

GEF=OSORT<GEF) 

DO 30 I = 1,NN 
DO 30 J = 1,2 


10 


15 


20 


d, J)>UH(K)<FGCOS( I, J,K) 
(dJI^-UHlKI^ACOSI I,J,K) 


32 





uuu 


t 


4 


35 

40 


45 

50 


55 


DO 3C K = 1,3 

DO 40 J = 

ABH(I,J| » 0.0 
00 35 K = 1,3 

ABHCIvJ) = ABH(I,J) ♦ HVP(I,J,K)4*2 

^5”ii*J»=DS0RT(ABHCI,jn '**''*'"* 

00 50 I — 1,NN 

cEin = 0.0 

00 45 J = 1,3 

CE(I> = CE(I) ♦ HVP<1,1,J)*HVP(I,2,JI 

NN2 i 

NNTT = 2*«{24NN| 

HNNT = 244NN 
M 55 I J l,NN2 

IC = 0 

DO 70 I = 1,NNTT 
IC = IC*1 
DE(IC) = 0.0 
FINTdC* * 1,0 
00 60 J = 1,NN 
J1 == IQ(J| 

J2 = lO(NN^J) 

oilic! : of»DE(IC?''^‘' * StJ2l«ABHIJ,2| 

FINTJICI = FINTIICI/HNNT 

IfTiO(J»-2» 66,65,66 

J * J+1 
GO TO 64 
IQ(J) = 2 
CONTINUE 

FORMAT (IH ,10F11,5) 

DO 200 I*1,IC 

* De(n/(1.39963*GEFI 

END 


60 


64 

65 


66 

70 

100 

200 
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o 
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subroutine TWSF(A,B,ACOS,OECOS»GCOS»GECOSfPAR» 

THIS PROGRAM SETS UP DIRECTION COSINES FOR THE PLANAR 
RADICAL IN THE CRYSTAL AXIS SYSTEM. THEN THE ANGLES FOR 
TWIST, WAG, FLAPf AND HCH ANGLE ARE USED TO RE-ORIENT THE 
RADICAL. THE DC*S ARE CHANGED FOR A CHANGE IN ANY OF THESE 
PARAMETERS. 

Ad) AND Bd) - KNOWN VECTORS IN THE ORIENTED RADICAL. 

IMPLICIT REALMS<A-H,0-Z) 

DIMENSION M3).B<3),YO).CHlf3),CH2I3I.XlC3),X2(3),ACOS(2,3,3) 

HHCH=PAR(4)/2.0 
CANG=90.0-HHCH 
HHCH=HHCH*3.1416/180.0 
CANGsCANG43.1416/180.0 
D»OCOSfHHCH) 

E^^OCOSICANG) 

DO 10 1=1,3 
CHin)=04A( Id-E«Bd ) 

10 CH2U)=D*A(I)-E»B(I ) 

CALL CR0SSP(Y,CH1,X1I 
CALL CR0SSPfYtCH2»X2) 

DO 15 1=1,2 
DO 15 J=l,3 
DO 15 K=l,3 
6EC0Sd,J,K)=0.0 
15 GC0S(I,J,K)=0.0 
J=1 
K=2 
L=3 

DO 20 M=l,2 
DO 20 1=1,3 
GO TO f30,40),M 
30 ACOS(M,J,n=CHlU) 

ACOS(M,K,I)=Yd) 

AC0S(M,L,I)=X1U) 

6C0S<M,J,I)=A(n 

GCOS(M,K,I)=B(!) 

GCOS(M,L,n=Yd) 

GO TO 20 

40 AC0STm,J,I)=CH2(I) 

ACOSlM.K,I)=Yd) 


r 


4 
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< 


-h 


20 


50 


ACaS(M,Lf n==X2U) 

CONTINUE 

CALL AXR(A,ACOS.PARCin 
CALL AXR(Y,AC0S,PARC2») 
CALL AXR(B»AC0S*PARf3)l 
CALL AXR(A»GC0S»PAR(1)) 
CALL AXR(Y,GC0S»PAR(2)I 
CALL AXR(BtGC0StPARf3n 
DO 50 1=1,2 
3 
3 

K)=AC0SfI,J,K) 
K)=GC0Sf 


00 50 J=l, 

DO 50 K=l, 

DEC0S(I,J, 

GECOSd.J, 

Y(1I=0.0 

Y(2»=1.0 

Y(3)=0.0 

CALL AXRCY 

CALL AXR(Y 

RETURN 

END 


,GECOS,180«) 
tOECOSflBO.) 


SUBROUTINE CROSSPfA,B,C) 

THIS PROGRAM CALCS THE CROSS PRODUCT BETWEEN TWO VECTORS* 

IMPLICIT REAL*8(A-H,0-Z» 

DIMENSION A(3»,B<3»,C(3» 

Cf 1)=AC2)4‘B(3}-A(3)«B(2) 

Cf2l=Af3}*B(1)-A(1I«B(3) 

Cf3)=A(l)*B(2I-'A(2l*B(l) 

RETURN 

END 


SUBROUTINE AXR(A,ACOS,ANGLE) 

PROGRAM CHANGES THE DIRECTION COSINES FOR EACH TENSOR FOR A 
CHANGE IN ANY PARAMETER f ANGLE) •‘uk CAcn i twiUK l-UR A 

IMPLICIT REAL*8CA-H,0-Z) 





ooo 


DIMENSION A(3»,ACOS<2,3,3»,ROTM(3,3),ROTMT(3,3l,AR(3,3I, 


o 


DO 10 1=1,3 
DO 10 J=l,3 
10 ROTM(I,JJ=0.0 
R0TM(1,1»=A<2I/E 
R0TM(1,2)»-A(1»/E 
R0TM<2,1)=A(1»*A(3)/E 
R0TM(2,2I=A<2»*A(3I/E 
R0TM(2,3»=-E 
R0TM(3,1)=A(1) 

R0TMC3,2)=A<21 

R0TM(3,3)=A(3| 

CALL GMTRA(R0TM,R0TMT.3.3» 
ANGN=ANGLE*3«14159265/160. 
F=0C0S(ANGN) 

G=0SIN(ANGN) 

00 20 1=1,3 
DO 20 J=l,3 
20 AR(I,J)=0.0 
AR(1,1)=F 
AR(l,2i=G 
AR(2,1)=-G 
AR(2,2I=F 
AR(3,3)=1.0 

CALL GMPRD<AR,ROTM,S.3,3,3) 

CALL GMPRD<R0TMT,S,Sl,3,3,3) 

L=2 

DO 30 J=l,3 
DO 30 K=l,3 
OCOSKK, J)=AC0S(I ,J,K) 

30 0C0S2<K,J)=AC0S(L,J,K) 

CALL GMPR0(S1,0C0S1,DC0S3,3,3,3) 
CALL GMPR0IS1,0C0S2,0C0S4,3,3,3) 
DO 40 J=l,3 
DO 40 K=l,3 

AC0Sn,J,K)=0C0S3<K,J) 

40 ACOSIL,J,K):b0C0S4<K,J) 

RETURN 

END 


« 




S(3,3I 











105 


C 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


A GENERAL MATRIX. 


SUBROUTINE GMTRA(A,R,N,M| 

THIS PROGRAM WILL TRANSPOSE 

IMPLICIT REAL*8(A-H,0-ZI 
DIMENSION A<1)«R(1I 
IR=0 

DO 10 I-1,N 
IJ=I-N 
DO 10 J=lfM 
IJ=IJ+N 
IR=IR+1 

10 RCIR)=AfIJ) 

RETURN 

END 


SUBROUTINE GMPRDIA»BtR,N,NtL) 

THIS PROGRAM WILL CALC THE PRODUCT OF 2 MATRICES, 

IMPLICIT REAL*8(A-H,0-ZI 
DIMENSION AU>tB(l)»R(l) 

IR=0 

IK=-M 

DO 10 K*1,L 

IK=IK+H 

DO 10 J-1,N 

IR*IR+l 

JI=J-N 

IBsIK 

RCIR)=0 

DO 10 1=1,M 

JI=JI+N 

IB=IB+1 

10 R(IRI=R(IRKA(JII«Bf IB) 

RETURN 

END 


A SAMPLE OF THE DATA IS AS FOLLOWS 
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.575092 

•269381 

0.0 

-61. 
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1 2 
1 2 


2 13.70 
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.455907 

.579567 

0.0 

-91. 

3 12 

3 12 

3 12 


1 . 1 . 
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0.0 

2.0042 
2 3 1 

2 3 1 

2 3 1 
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2 

2 

2 


3 

3 

3 
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1 
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.70711 0.0 
1 16 
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1 

1 

1 

0 

0 
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3 

1 

1 

1 

0 
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1 

1 

0 
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3 

1 

1 

1 

0 
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3 

1 

1 

1 

0 
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1 

2 

1 

0 
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3 

1 

2 

1 

0 
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3 

1 

2 

1 

0 
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• .0001.0 
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- 
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